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LINEAR ALGEBRA USING PASCAL MT+ 

Introduction to Linear Algebra Calculator 

This program solves the two major problems of linear 

algebra: (1) the determination of the eigenvalues and 

eigenvectors of a square matrix "A" and (2) the solution of 

simultaneous equations "A X = B" for an n by n invertible matrix 

"A" and any n by m matrix "B" . One must use the built-in matrix 
editor to enter the matrices "A" and "B" . A program , based on 
eigenvalues of a companion matrix, to calculate the roots of a 
polynomial with real coefficients is also included. The source 

code of the mathematical routines and I/O routines could be used 

to build a more elaborate matrix calculator that might include 
disk l/O , a more flexible editor, addition and multiplication cf 

matrices, and a matrix expression evaluator. This program does 

not have these capabilities. 

The hardware configuration necessary to run this program is 
an IBM PC with 128I< of RAM memory and one double-sided drive. A 
is written in Pascal MT+ and thus the source code could be 
complied with any MT+ compiler after making a few modifications 
to the system I/O routines. Depending upon what floating point 
library is linked in with the compiled code, the program may be 
run with an 8037 chip or an 3088 chip. 

To run the program one may do any one of the following: (1) 

if the machine is off, insert the progam diskette in arive A _nd 
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turn the power on; CPM-86 will be loaded and the program, called 
MATRIX, will be executed; (2) if the machine is on, insert the 
diskette in drive A and do a warm boot by pressing "Ctrl-Alt-Del" 
simultaneously; CPM-86 will be loaded and the program will be 
executed; or (3) if CPM-86 has been loaded, insert the diskette 
in drive A, make the CPM prompt say A>, and key in MATRIX 
followed by carriage return. There are obviously other 
possibilities. The program and all its overlays could be 
transferred to RAM disk M and executed from there with M also the 
default drive. One must always have the program and its overlays 
on the default drive because the program looks for its overlays 
on the default drive. 

Brief instructions on how to use the program will now be 
given. When the program is executed, the main menu is displayed. 
By pressing "G", one can view some very brief instructions. They 
are intended for somebody who has never run the program before or 
who does not like to read documentation. Before doing any matrix 
work, one must press "F" to bring up the matrix editor ana 
displayer. The menus are self-explanatory and channel the user in 
the desired directions. To do matrix work one must enter a matrix 
"A" but tne entering of "B" is only necessary for solving "AX = 
B" . Only one matrix "A" and one matrix "B" may be stored at a 
time. To enter a new matrix after one has already been created 
one must answer "Y M to the erase prompt. 

When matrices "A" and "B" are first created, they are 
initialized at all zeros. The user then edits the matrix in 
increasing column order. Each column will be displayed for 
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editing. If one answers "Y" to the column correct prompt, then 
the next column will be displayed for editing. If one answers 
"N", the cursor will be positioned at the entry in the first row 
of that column. If one wants to change that entry, he enters u 
now number followed by return. If one simply enters return, the 
entry is not changed and the cursor is positioned at the next 
row. After the entire column has been edited, the column is 
redisplayed with the same prompt. To edit column 4, one is forced 
to first view and answer "Y" to the first throe column correct 
prompts. A compromise on flexibility and editing speed for user 
friendliness and less programming has been made here. 

Once "A" has been entered, one may return to the outer menu 
to calculate the inverse, determinant, eigenvalues, and 
eigenvectors. If "B" has been entered also, one can press "C" to 
solve a linear system "A X = B". One can solve for the roots of u. 
polynomial by pressing "A". One enters the coefficients of the 
polynomial within this segment of the program, not in the 
editor . 



Two possible I/O problems are the following: (1) real numbers 
are not filtered for underflow or overflow when inputted; they 
have to be entered in the interval [1*0 E-3U7, 1.0 11+307J; and 
(2) the printer must be on-line ( or must be put on-line when the 
CPrt-86 error message is displayed) if a hardcopy is desired. If 
one cannot get the printer on line after a hardcopy was 
requested, CPM-86 exits the user back to a system prompt and the 
entered matrices "A" and "B" are lost. 
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An example matrix "A" to key xn is : column 1 is the vector 
(1,3,6), column 2 is the vector (-3, -5, -6), and column 3 is the 
vector (3,3,4). Let X and r (not variables in the Program Matrix) 
stand for the associated eigenvector and eigenvalue in the 
equation "AX=rX" . The eigenvalues of "A" are 4,-2, -2 and the 
associated eigenvalues are respectively (1,1,2), (-1,1,2) and 
(1,3,2). It should be noted that eigenvectors returned by the 
calculator will differ by scalar multiples from those above. 
Also, keep in mind instead of (-1,1,2) and (1,3,2) which are 
associated with the eigenvalue -2, we could have -another pair of 
eigenvectors such as (0,4,4) and (2,2,0) which are in the same 
eigenspace as the first pair. The inverse of the matrix "A" is 
the matrix whose columns are: (1) . 125 ( -1 , 3 , 6 ) , (2) .125(-3,-7,- 
6), and (3) .125(3,3,2). This matrix has 16 as its determinant. 
The polynomial X**s-oX**2+llX-6 would be entered as the vector 
( 1 , -6 , 1 1 , -6 ) . It has roots 1,2,3. 

Description of CRT and keyboard utilities 

The utility programs in Module CRTLIB are for cursor 
control and filtered input of characters, integers, and real 
numbers. This section should be read concurrently with the source 
listing of Module CRTLIB. These routines are based upon two CPM- 
86 BIO J . ills, which are implemented in the Procedure Bioscall. 
Bioscall has two input parameters: (1) FUNC, a variable of type 
CPHOPERATION and (2) 0CI1, an integer variable. Bioscall, as 
presented nere, only responds to two inputs, COLn'IN and CONUUT, of 
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type CPMOPERATION. 



In the following, terminal character refers to a character 
that is being sent to the monitor to effect a cursor control. To 
send a terminal character to the CRT interface ( monitor ) , one 
calls Bioscall with input parameters FUNC equal to CONOUT and OCH 
equal to the ASCII number of the terminal character. The value of 
the variable BDOSVAL( integer global to Module CRTLIB) is not used 
explicitly; here it is used only to call Procedure 0BDOS86. To 
read a character from the keyboard, one calls Bioscall with input 
parameters FUNC equal to CONIN and OCH equal . to any dummy 
integer. The global variable BDOSVAL is set equal to the ASCII 
number of the keyboard character pressed. The number 50 in 
@BDOS86(50, ADDR ( DESCRIPT L 1 J ) ) is the BDOS function number for a 
BIOS call. The five consecutive bytes allocated by DESCRIPT are 
used for passing information in the 9BDOS86 call. For example, 
"DESCRIPTl 1 j : =4 " is the Pascal line for console display while 
"DESCRIPT[1] :=3" is the line for keyboard input. The Pascal MT+ 
manual gives a brief description of the MT+ utility (§BD0386 while 
the CPM-3G manual gives a detailed discussion of BDoS calls in 
its Appendix D. The Procedure Bioscall presented here is a Pascal 
implementation combining 0BDOS86 with the BDOS calls (for function 
number 50) of CPU-86. 

One could expand Procedure Bioscall to include all the 
parameters in CPMOPERATION by combining the ideas in Bioscall 
with Appendix D of the CPM-86 manual. One does not need to use 
Bioscall to write an alpha-numeric character to the screen; Write 
or Wnteln can be used for this. Most of the routines in Module 
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CRTLIB either call Bioscall or another routine which calls 



Bioscall. Function Getchar is the only procedure that calls 
Bioscall with the CONIN parameter. 

in order to use the utilities in Module CRTLIB one must call 
Procedure Crtinit at the beginning of the main program. Crtinit 
initializes the arrays CRTINFO and PREFIXED so that their values 
can be used by the utilities in Module CRTLIB. 

Function Getchar, using Bioscall, performs the task of 
reading a character from the keyboard. The input to Getchar is a 
variable set of characters, called OKSET in Getchar 1 s 
declaration. When Bioscall is called by Getchar to get a 
character from the keyboard, Getchar then checks to see if the 
character represented by BDOSVAL is in OKSET. If it is, this 
ciiaracter becomes the value of the Function Getchar and is used 
in the procedure calling Getchar; if it is not, a beep is sounded 
and the process is repeated until a character in OKSET is finally 
entered. This version of Getchar automatically changes lower case 
to upper case. An example of Getchar being used with OKSETs that 
are dynamically changing can be seen in Procedure Getreal. 

Procedure Crt is tne main routine that calls Procedure 
Bioscall with the CONuUT parameter. Depending on the boolean 
value given by an element in the array PREFIXED, the leuum 
character, ESC, is sent to the monitor and then the actual 
terminal ciiaracter is sent to the monitor. For example, suppose 
we want to clear the screen. One uses the variable ERASE03 of 
type C RT CO OMAN D as the parameter input to Procedure Crt. Using 
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ERASEOS as the index variable, Crt sends the ieadin character 
ESC to the monitor since PREFIXED[ERASEOS] is true and then sends 
the character ' J', which is CRTINFO[ERASEOS ] . It should be noted 
that in the above, the ASCII numbers of the various characters 
are actually inputted to Bioscall in Procedure Crt. 

One can control background colors by using the type COLOR, 
the arrays COLORIHFO, CRTINFO, and PREFIXED along with the 
Procedures Initcolor, Altcolor, and Paint. One then needs to 
write a procedure similar to Altcolor to change foreground 
colors. These would be the fundamental library routines for using 
colors . 

Procedure Gotoxy (X, Y: INTEGER ) , which also calls Bioscall 
with input parameter CONOUT, places the cursor at vertical line 
number X and horizontal line number Y. The other cursor control 
routines are fairly obvious in view of the above discussion. 

Procedure Intread is used to read an integer between -32768 
and +327o7. The characters are filtered, put into a string, 
checked for the proper range, and then converted to an integer. 
Intread requires an integer to be entered; one cannot simply 
enter carriage return. Procedures Getreal and Value are used 
together to read a real number. Some explanation is required on 
the use of Getreal. In the procedure calling Getreal, we have a 
string variable, SREEL, initialized at ' ' , an empty string of 
length zero but still a string. Procedure Getreal reads and 
dynamically filters a string of characters and if this string is 
of length greater than zero, this string is returned to the 
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calling procedure by Getreal and SREEL is set equal to it. If 



Getreal only reads a carriage return, then SREEL remains only an 
empty string and Procedure Value is not called. If SREEL lias 
length greater than zero. Procedure Value calculates the real 
number corresponding to SREEL. Getreal and Value are used to edit 
a real number that already exists, either through initialization 
or calculation. By pressing carriage return, one can leave an 
existing real number as is. This use eliminates a computer crash 
caused by keying in only a carriage return when the usual READ or 
READLN wants an actual number. Procedure Value was- taken from The 
Byte Book of Pascal while Procedure Powrten is from the Pascal 
Users Manual by Jensen and Wirth. 



Data and Overlay Description 

This program illustrates the use of module overlays, an 
ordinary module, matrix types and external variables. This 
description explains the manner in which the program was 
implemented, and therefore it should oe used as a guide to 
understanding the structure of the code. Modules are separately 
compiled groups of procedures. There are two types: 11) an 
overlay module and (2) an ordinary module which must be linked 
with an overlay module or main program. The Module CRTLIB is 
separately compiled and linked in with the main program. CRTLIB' s 
subroutines essentially become a part of the main program! 
sometimes called the root overlay) and always are resident m RAM 
memory. Since the utilities in CRTLIB are continually being 
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called, it would not be practical to make CRTLIB an overlay 
module and thus force a lot of disk accesses. If a module 
overlay ( including the main program) uses a procedure in Module 
CRTLIB, then the "calling" overlay must declare the procedure as 
an external procedure in the calling overlay's declaration 
heading just below the variable declarations. 

The data storage of the program will now be explained. The 
controlling data allocation declarations are (1) the TYPE 
DOMAIN 1 = 1 . . 20 , (2) the CONST MAXDEG=2G, (3) the CONST 
MAXDEGP1=21, and (4) the TYPE D0MRP1=1 . .MAXDEGP1 . MAXDEG, 
MAXDEGPi (=MAXDEG PLUS 1), and DOt-lRPl are used to solve for roots 
of a polynomial. D0MAIN1 is the maximum dimension of a matrix 
entered as input data. MAXDEG should De set to be less than or 
equal to the value of the D0MAIN1 upper bound since the 
polynomial root finding technique uses an n by n matrix to solve 
an n degree polynomial equation. In this program, the 
coefficients of a polynomial are not stored in a global variable! 
for the main program) and there is not much I/O to go with the 
polynomial routine. It is included as a side benefit to 
illustrate the use of the eigenvalue procedures as subroutines to 
another procedure. Variables that are referred to as global are 
declared in the main program and then declared external in the 
overlays that use them. The matrix "A" whose square dimension 
is IOMDIM is stored in the upper left-hand corner of the global 
variable IOTA and the matrix "B" whose dimension is ROWDIMB by 
IOM is stored mthe upper left-hand corner of the global variable 
I0TB. IOTA, 10TB , I0NDIM, R0WDIM3, and IOM are declared m the 
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root overlay! Program Matrix) and referenced in the appropriate 
overlays by declaring them as external variables. The global 
matrix "XX" is used to store the solution to "AX=B" and is 
referenced as external in Module Overlay2U( whose source code is 
AXEOSB.SRC). The global vectors TEVR, TEVI store the real and 
imaginary components of the eigenvalues and are referenced as 
external in Module Overlay24( EIGENVEC . SRC ) . The global matrix 
TVEC, referenced as external in Module Overlay24, is used to 
store the eigenvectors of the matrix "A" . Communication of data 
between the root overlay (Program Matrix) and the overlay 
modules if 3 , # 19, 42V , ‘424 is accomplished by initially declaring the 
above mentioned variables in the main program and by declaring 
them external in the necessary overlay modules. The Types of 
matrices and vectors must be declared in the root overlay and 
redeclared in the overlay modules. When using dimensions less 
than the maximum allocated in the Types, all data will be stored 
m the upper left-hand corner of matrices or beginning components 
of a vector. Intermediate results produced by Module Overlay 1 
( LIMSYS . SRC ) , Module Overlay5 (EIGENHQR. SRC ) , and Module Overlays 
t, EIGENBAL . SRC ) are passed as parameters in the procedure calls. 

The root overlay (Program Matrix) ea Lis airectiy six main 
procedures located in six overlays. They are (1) Procedure 
Matrixio, located in Module 0verlay3, which is the editor for 
matrix input, (2) Procedure Help, located in Module Overlay23, 
which is a brief set of instructions on how to run the program, 
(3) Procedure Ttrdet, located in Module 0verlayl9, which 
calculates the determinant of the matrix stored in variable IOTA, 



(4) Procedure Ttsaxb, located in Module Overlay2d, which solves 
"AX=B" for "A" stored in the variable IOTA and "B" stored in the 
variable IOTB or "B" is the identity matrix ( used to calculate 
the inverse of "A"), (5) Procedure Ttrnaa, located in Module 
Overlay24, which calculates eigenvectors and eigenvalues of the 
matrix stored in IOTA, and (b) Procedure Ttrpqr, located in 
Module 0verlay21, which calculates the roots of a polynomial with 
real coefficients. These six procedures which are located in 
module overlays must be declared external in the root overlay 
heading. After the word "External", the overlay number to which 
the procedure belongs must appear in square brackets. 

Other external procedures used witnin these overlays do not. 
have to be declared within the root overlay; they are declared as 
external procedures with the appropriate overlay numbers within 
the overlay that calls them. For example, Procedure Ttrdet, 
located in Overlay Modulel9, calls Procedure Rlud, located in 
Overlay Modulel . Therefore, ' Procedure Rlud must be declared 
external in the heading of Module 0verlayl9. Module Over lay 3 uses 
procedures located in the Module CRTLIB( not an overlay module), 
ana therefore these called procedures must be declared as 
external m the Heading of Module Overlay3. These external 
procedures from Module CRTL1B do not have overlay numbers. 

Module Overlay3 ( IOMATRIX . SRC ) is the module responsible for 
keyboard input of matrices and printer output of matrices. The 
maximum size of the matrix is set m the TYPE DOMAIN 1 = 1 . . 2b . If 
one wanted only matrices of dimension lb or less, one would set 
TYPE DOMAIN! = 1 . . lb . One must change these type declarations for 
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all the overlays that use them. It should be noted that che 
editing programs are based upon editing a column of a matrix that 
fits onto one screen display. If one wanted to enter a column of 
dimension larger than twenty, the column display and editing 
routines would have to be modified a little. On the other hand it 
is not practical to enter a matrix larger than 10 by 10, let 
alone 2u by 20. Large matrices should be generated by a program 
or read from a disk file. For example. Module Polyroot generates 
a 20 by 20 matrix when one requests the root of X**20=l. 

Within this module, Procedure Getcolumn and Procedure 
Getpolcof are the work horses. Respectively, they take the column 

y 

of a matrix or the coefficients of a polynomial and display them 
on the screen for editing. These two routines use the cursor 
control utilities and Procedures Intread, Getreal, and Value to 
edit existing columns. Procedure Matnxio is the editing 
procedure which calls the Procedures Edit and Display. Edit 
displays menus, initializes matrices, and calls the above 
mentioned procedures to edit existing matrices. Procedure 
Display, which is somewhat redundant, gives a faster display than 
Procedure Edit. Procedure Hardcopy is used to printout the 
matrices "A", "B", inverses, determinants, matrix equation 
solutions, eigenvectors, and eigenvalues. 

Module Uverlay24 (EIGENVEC . SRC ) , Module OverlayG 
( EIGENBAL . SRC ) , and Module 0verlay5 v EIGENHQR . SRC ) are the 
overlays which calculate eigenvectors and their corresponding 
eigenvalues. Procedure Ttrnaa takes the global matrix IOTA 
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containing the matrix "A" and its dimension, IONDIM, and calls 
the Procedure Rnaa . Procedure Rnaa is the procedure which calls 
the Procedures Balance, Elmhes, Eltran, Hqr2, and Balbak, which 
together calculate the eigenvalues and eigenvectors. The 
eigenvalues are returned to Ttrnaa in the global arrays TEVR and 
TEVI . TEVR[ I j and TEVi[l] give the real and imaginary parts of 
the eigenvalue #1. The eigenvectors are returned in the matrix 
TVEC in a not so obvious way. They are returned in an order 
corresponding to the eigenvalues in TEVR and TEVI. The manner in 
which TVEC returns the eigenvectors needs to be illustrated by an 
example. Let TEVRL1J=1 and TEVIL1 j= 2; the eigenvalue is l + 2i. 
Then TEVR[2]=1 and TEVl[2j=-2 since eigenvalues and eigenvectors 
of real matrices occur in conjugate pairs. Now the two 
eigenvectors that correspond to these eigenvalues are conjugates. 
TVEC in its first column will have the real part of eigenvector 
#1 and in its second column it will have the imaginary part of 
eigenvector #1. Eigenvector #2 will just be the conjugate of this 
eigenvector. If the eigenvalue is pure real, then the eigenvector 
is pure real and only one column is necessary. In the case of a 
matrix with repeated eigenvalues without a full set of 
eigenvectors, the eigenvectors will be repeated in the matrix 
TVEC( probably numerically slightly different). These eigenvector 
subroutines are Pascal translations by the author of LIdPACK 
eigenvector routines in Fortran. 

Module Overlay21 l POLY ROOT . SRC ) also uses the eigenvalue 
subroutines to solve for the roots of a polynomial with real 
coefficients. It makes the integer 1 the leading coefficient of a 
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polynomial of degree n and enters the other n coefficients into 
the first row of an n by n matrix with l's on the subdiagonal; 
this matrix is called the companion matrix and the eigenvalues of 
it are the roots of the polynomial. Procedure Ttrpqr does the 
preliminaries and Procedure Rpqr calls the sequence of eigenvalue 
subroutines to find the eigenvalues wnicn are returned in the 
variables TWR and TWI . Procedure Balbak is not used since the 
eigenvectors are not needed. POLYROOT. SRC illustrates the use of 
these subroutines in a real problem and it also serves as a test 
of the correctness of the eigenvalue codes. 

Module 0verlay2d (AXEQSB.SRC) and Module Overlayi 
(LIN3YS.SRC) are the overlays which solve systems "AX=B" as well 
as find the inverse of the matrix "A" (B will be the identity 
matrix m this case). Again the "A" matrix is in IOTA and the "B" 
matrix is in IOTB. Procedure RX.UD does an L-U decomposition on 
the matrix in the variable SA( the constant matrix SB tags along) 
and Procedure Rfbs does the back substitution to solve the 
system. In the case where "B" has more than one column, Rfbs is 
called IOM times to solve the equations. Procedures Rlud and Rfbs 
are Pascal translations by the author of a Sandia Labs matrix 
equation solver; it is very successful with matrices "A" that are 
almost singular(no inverse). 

Module 0verlayl9 (DETER11. SRC ) and Module overlayi 
(LIN3YS.SRC) are used to calculate determinants. Rlud in 
LINSYS.SRC does an L-U decomposition of the matrix TALU ( a copy 
of IOTA) and returns the L-U decomposition in TALU. The matrix L 
has ones on the diagonal and therefore the determinant of "A" 
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will be the determinant of "U" which is just the product of U's 
diagonal entries. 



Layout of Overlays 

Module overlays are divided into groups. Overlays ffl -ffJLG 
are in overlay group #1, overlays #17-#32 are in overlay group 
#2, etc.. The source code of an overlay module must begin with 
its identif ication number; for example overlay module #1 t,. 3 the 
heading, Module Overlayl . The root overlay is just the main 
program; here it is Program Matrix. This program does not use the 
most sophisticated overlay structure. In this program, overlays 
within the same group do not call each other and there is no 
heap. When compiling an overlay, one does not have to specify the 
address in RAM memory where the overlay will be loaded. This is 
done during the linking process. All the overlays within the same 
group must be located at the same address. This address must be 
larger than the length of the root overlay including any runtime 
library routines or ordinary modules linked with tne root 
overlay. For example, the root overlay, Program Matrix with 
routines from the libraries TRANGEND . R86 , 87 REALS . Rob , 
PASLIB.R86, and Module CRTLIB use bFFCH(hex) bytes of storage. 
Tile memory allocated for overlay group £ 1 must be some address 
larger than this; here 36UUH is used. In tiie linking instructions 
below, overlays £1 - £16 must be directed to begin at 8600II. Then 
one looks at the largest length of the compiled overlays in group 
#1 and adds this to 860011. In tills program, Module Overlays has 
length 42AAH. This should be verified again after all the 
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overlays in group #1 have been linked to the root overlay. One 
adds 8600H to 42AAH to obtain the length C8AAH. The address where 
overlays £17 - £32 begin must be larger than C8AAII. This program 
uses address D200H for the location of overlays in group if 2. Tne 
overlay of the largest length in group £2 is Module Overlay21 
with a length of 1370H. 137QH plus D200H equals E570H. Therefore, 
the total length of the code that must be specified in the last 
link command with the R switch must be greater than E570H. The 
value used is FC0GH. The data storage required is only that from 
those variables declared in the main program and it is 57GCH. 
The value 8U00H is more than adequate and is declared with the D 
switch in the last linker command. Local variables and parameters 
are stored in the stack. This stack is automatically given 32K 
bytes of RAM. One can actually assign the stack size by using the 
Z switch in the last linker command. For example, "Z/3u b" wili 
assign 30U0H for stack memory. Note that 300 here means 3bOUH or 
12K bytes. 

One has to make sure the total RAM allocated for data, code, 
stack, and heap is smaller than the RAM allocated for program 
execution. Otherwise, upon program execution, a "Memory not 
available" message will be displayed. In this program, the stack 
was reduced to 12K in order that the program would execute in a 
128K machine. Having a stack too small is a potential problem, 
but by far the major source of problems is having one overlay 
group overwrite another. Apparently, one needs to have some 
memory available between the end of the largest overlay in a 
particular group (including the root overlay) and the beginning 
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of the next overlay group. For example, in this program when the 
addresses 7200H, B800H, and DE00H were used in place of the 
addresses 8600II, D200II, and FC0UH the program would not execute 
properly. The Pascal MT+ manual should address this problem. If 
one experiences weird program execution, the distance between 
overlay groups should be increased as the first approach to solve 
the problem. 



Computer Configuration for Pascal MT+ 

in doing developmental work for this program, the IBM PC was 
configured with two DS/DD disk drives, 512K of RAM , and the CPM- 
86 operating system. The RAM was partitioned into 14UK for 
running programs ( compi ler , editor, linker, and the executable 
linear algebra routine) and J72K for a RAM disk, called the M: 
drive. The program SETUP.CMD, furnished with CPM-86, was used to 
configure the RAM drive and the cursor pad of direction arrows. 
Upon bootup, the RAM drive is created and the cursor pad is 
initialized . 

A RAM disk is essential (though not required) for program 
development with Pascal MT+ because of the heavy disk access by 
the MT+ Speed Editor, Compiler, and Linker of their overlays and 
the program's source and object code files. If one follows the 
implementation techniques described in the following, the user 
will find MT+ not just a powerful Pascal implementation, but one 
that possess good development time and is not cumbersome. It 
should be mentioned that in many articles comparing the MT-r 
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implementation to others, the benchmarks involving compile and 
link times are misleading. Pascal MT+ uses a lot of disk I/O 
because of its code size and this eats up a lot of time. Using a 
RAM disk reduces this I/O time considerably. 

One must configure three diskettes. Diskette frl, called the 
compiler diskette, must contain CPM.SYS( configured for RAM drive 
and optionally the cursor pad), the MT+ compiler ( consisting of a 
CMD file and its several overlays), and three R86 files: (1) 
8 7 REALS . R8G , (2) PASLIB . R86 , and (3) TRANCEND . R86 . It is optional, 
but useful to put the utility routines such as PIP. CMD and 
STAT.CMD on the compiler diskette. Diskette if 2 , called the editor 
diskette, should contain the Speed Prograiamming Editor ( CMD file 
and its overlays;, SIP.CMD( PIP. CMD renamed by the user), ana 
STAT.CMD. Diskette #3, called the linker diskette, should contain 
the MT+ I inker ( CMD file and its overlays) and the SUBMIT.CMD 
file. For single stroke execution, it is convenient to rename the 
compiler, editor, and linker respectively M.CMD, S.CMD, and 
L . CMD . 

The following is the start up’ process. Put the compiler 
diskette in drive A: and the editor diskette in drive B: and then 
turn the power on. When the "A>" prompt comes up, key in "PIP 
M:=B:S*.*". This will transfer the editor, SIP. CMD, and STAT.CMD 
to the M: drive. One can use an editor other than DRI ' s Speed 
Editor. Users of UCSD's P-Editor will find the Speed Editor 
acceptable when not in the P-system. After the editor is 
transferred to RAM drive, one removes the editor diskette from 



drive B: 



and inserts a work diskette containing source and their 



compiled, R86, files* 



One then presses Control-C and changes the default drive to 
Ms . if one wants to edit an existing source file, one uses 
SIP.CMD on drive Ms to "pip" the source file from drive 3: to 
drive Ms. While in the Ms drive, one edits a Pascal source file. 
If one is using the Speed Editor, a word of caution is necessary. 
If one attempts to save a file on a drive without enough memory ( 
including RAM drive), the old source file as well as the source 
file in the editor's buffer will be lost. This potentially 
serious problem can be avoided if one is aware of it. If the 
Speed Editor is used one should use its syntax scanner and 
variable checker to locate syntax errors and typos prior to 
compilation. 

Once a source file has been edited it must be compiled. It 
is not necessary to transfer the compiler to RAM drive. One must 
change the default drive from Ms to As, leaving the compiler 
diskette in drive A: and the source file on RAM drive. Since 
Pascal MT+ is oriented toward modularity and modular compilation 
in particular, individual source files( modules and overlay 
modules ) are relatively short. Because of this, disk I/O to 
drive A: to load the compiler overlays is not too frequent. With 
the source code being read from RAM drive and compiled code being 
written to RAM drive, compiles are quite timely. 

To compile all the overlays and modules for the linear 
algebra calculator, one must do the following tasks. There are 
eleven source codes which must be saved on the work disk. They 
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are (1) MATRIX. SRC, (2) CRTLIB.SRC, (3) IOMATRIX . SRC , (4) 
EIGENBAL . SRC , (5) EIGENHQR . SRC , (6) MHELP.SRC, (7) LINSYS.SRC, 
(8) DETERM. SRC, (9) AXEQSB.SRC, (10) POLYROOT . SRC , AND (11) 
EIGENVEC . SRC . Each file should be transferred to RAM drive and 
then compiled. To compile the root overlay, MATRIX. SRC, one 
enters the following command: "MT+86 M: MATRIX . SRC" . If one 
changed the name of the compiler to M.CMD, then one enters "M 
M : MATRIX. SRC” . The compiler will compile the file and write a 
code file, MATRIX. R36, to RAM drive. Similarly, one compiles the 
other ten source code files to produce ten more RSb coue files. 
These code files should be saved on the work diskette in drive 
13:. It should be emphasized that if the source code of a module 
or overlay module is not changed, then it does not have to be 
compiled . 

The most important implementation method in MT+ code 
development is to automate the linking process. Without this 
automation, the linking process is extremely time consuming. 
Automatic linking is accomplished by using SUBMIT.CMD and linker 
input command files, called KMD files. These will now be 
described in detail and applied to the linear algebra calculator. 
Once all the overlays have been compiled one makes up a KMD file 
for each compiled overlay module. A KMD file is actually a text 
file keyed in using an editor and given the suffix KMD instead of 
the usual SRC suffix. The following KMD text files are needed: 

(1) Md.KMD with the line of source: 



M : MAT RIX , M: CRTLIB , M : T RANGED D , M : 8 7 REALS , M : PAS LI B / S / E / W 



(2) Ml.KMD with the line of source: 



M: MATRIX. 00 l=M:MATRIX/0: 1 ,M:LINSYS,M: PASLIB/S/P: 8600/L 

(3) M3.KMD with the line of source: 

M: MATRIX. 00 3=M: MATRIX/0: 3 ,M: IOMATRIX,M: PASLIB/S/P: 0600/L 

(4) M5.KMD with the line of source: 

M: MATRIX. 00 5=M:MATRIX/0: 5, M : EIGENHQR, M : PASLIB/S/P: 8o0U/L 
(5J Mb.KMD with the line of source: 

M: MATRIX . 006=M : MATRIX/ O : 6 , M: EIGEN3AL, M : PAS LIB/ 3/P:8600/L 
(b) M23.KMD with the line of source: 

M : MATRIX . 017 =M : MATRIX/0 : 23 , M : MHELP , M : PASLI B/S / P : D2u0 / L 

(7) M19.KMD with the line of source: 

M : MATRIX . 01 3=M : MATRIX/0 : 1 9 , M : DETERM, M : PASLIB/S/P : D200/L 

(8) M20.KMD with the line of source: 

M: MATRIX . 01 4=M : MATRIX/ 0 : 20 , M : AXEQS3 , M : PASLIB/S/P: D200/L 

(9) M21.KMD with the line of source: 

M : MATRIX. 015=M: MATRIX/0: 21 , M : POLYROOT , M : PASLI B/S/P : D2uu /L 
vlO) M24.KMD with the line of source: 

M : MATRIX . 018=M : MATRIX/ 0 : 24 , M : EIGENVEC , M : PASLI B/S/P : D200/L 

2 i 



(11) MR.KMD with two lines of source: 



M: MATRIX, M : CRTLI B, M : TRANCEND, M : 87 REALS, /C 
M:PASLIB/S/D:8000/V1 :860D/V2 : D2U0/R: FCU0/Z : 300 

These linker input command files are stored on the linker 
diskette. As with the editor, RAM drive will be used extensively. 
One transfers the compiled, R86, files from the work diskette to 
the M: drive. Then one transfers the three R86 library files to 

tiie M: drive: (1) 87REALS.R86, (2) PASLIU.RSG, and (3) 

TRANCEND. R8G . If one does not use the 8087 chip then FPkEALS.RSO 
must be used in place of 87REALS . R86 . At this point, one should 
remove the compiler diskette, and insert the linker diskette. 
Then one enters Control-C and makes A: the default drive. The 

input command files save a lot of time since the overlay module 
names and addresses rarely change during program development. 
Remember, that the linker has been renamed L.CMD for keyboard 
simplicity. Instead of keying in: "L M:MATRIX,M:CRTLIB, M: 

TRANCEND, M: 87REALS, M: PASLIB/S/ E/W" one simply enters: "L Mo/F" . 

Without the input command files option, one would have to do 
eleven links, keying in by hand each time all of the names and 
addresses in the linker input command files. Instead, one merely 
executes the eleven commands: (1)"L MU/F", (2)"L Ml/F", (3)"L 

M3/F" , (4) "L M5/F", (5) "L Mo/F", (6) "L M23/F", (7) "L M19/F" , 

(8) "L M2 Li/ F" , (9) "L M21/F”, (10) "L M24/F", (11) "L MR/F" . 

one can further automate the above process by using 
SUBMIT . CMD to do batch processing of the above eleven commands. 
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To do this, one must prepare a textfiie, with the suffix, SUB. In 
the linear algebra calculator, this SUB file is called 
MLINK2.SUB. In this text file, each line of text will be a single 
command. The eleven command lines of this textfiie will be the 
command lines of the above paragraph. It should be noted that 
this SUBMIT . CMD code is a very sensitive code. When the 
MLINK2.SUB was constructed using the Speed Editor, the batch 
process would not execute. If MLINK2.SUB was constructed using 
WORDSTAR with no editing changes, it worked fine. Also, some 
editing such as exchanges gave no problem while deletes did. 
SUBMIT.CMD seems to perform better under CCPM-86 . 

The following is the file configuration to do an automatic 
link. All RBu files including compiled source files and run-time 
libraries must be on drive M: . The linker diskette, in drive 
A:, must contain L.CMD and its overlays, SUBMIT.CMD, MLINK2 . SJli , 
and the eleven linker input command files. 

With A: as the default drive, one enters "SUBMIT MLIWK2" and 
the linking process will be automatically executed. once 
automated, the linking process is not a big hassle. Linking MT+ 
code is the price one pays for modular compilation. The link will 
be fairly fast, considering the size of the codes. The fact that 
the R86 files are being read from RAM drive and the final 
overlays are being written to RAM drive creates a tremendous 
savings in time. The linker will create ten overlays: (1) 
MATRIX . CMD , (2) MATRIX. 001, (3) MATRIX. UU3, (4) MATRIX. UU5, (5) 
MATRIX . 006 , (6) MATRIX .017, (7) MATRIX . 01 3 , (8) MATkIX.ul4, (9) 



MATRIX. 01 5, 



(1U) MATRIX. U18. 



After tney have been created, they 



should be transferred to another diskette to be saved and 



executed as desired. It might be possible to put the linker and 
its overlays and the KMD files on RAM disk also if one had a 
larger RAM drive. It does not seem that there is a tremendous 
loss in time to load the linker overlays and read the command 
files from a diskette. If the R86 files are read from a diskette 
and the compiled overlays are wniten to a diskette, there is a 
tremendous loss of time. One has to keep in mind that SUBMIT.CMD 
is sensitive to the selection of drives that have the files that 
it is batch processing. It seems only experimentation tells what 
SUBMIT.CMD accepts; it went for the above configuration. 
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SOURCE LISTINGS 



(* VERSION C 23 5 *) 

PROGRAM MATRIX; 

CONST MAXDEG = 20; 

MAXDEGP1 = 21; 



YPE DCMRP1 = 1..MAXDEGP1; 

COMA I Ml = 1..20; 

MATRIX = ARRAY [DOMAIN 1, DO 
L I ST I = ARRAY [DCMAIN1] OF 
LISTR = ARRAY [DOMAINl] OF 
LISTRPl = ARRAY [DC 1RPI] 0 
CRTCOMMAMD = ( ERAS ECS , ERA 

FCOLOR, BCOL 
3LINKCN, 
3LIMK0FF) ; 
SET OF CHAF ; 



MAIN 
I NT 
REA 
F RE 
SEOL 
: R. , R 



1] 

EGE 



p D 

EVI 



OF 

R; 



, DC 
DON 



,a: 



'L ; 



. RIGHT, 
TV I DOFF 



LEFT , 
, IMTE 



LEAD 

NON, 



IN, 

T ; ]p 



TIME, 

E N OFF , 



SETGFCHAR = 

PTR = ''INTEGER; 

CPMCPERAT I ON = (CCLD3C0T, 

?UNCUT,RD 
DSKREAD , D 

STRINC4C = STRING [40]; 
STRINGS = STRING [u] ; 
STRINGSO = STRING [SO]; 



MAR M 
FIN, 
S K r . "R 



300 

HCM 

ITE 



- r~ 

e,c: 



N3TAT,cc::iN,cc:;ou r 

;ldsk ,SET. 



r\ , . 



VAR I , J , IONDI I, 10 I, ROMDIMD : INTEGER; 

f TL 7 CT * ^ 1 1 T' r • 

GKCET : 3ETCE CHAR; 

AC I PL AG , IOBFLAG , QUITFLAG , SY.MFLAG : 3 C C L L A : 7 ; 
ICT A , TA , TVEC : . 1 AT R I X ; 

ICT3 , TEMPTS , XX : MATRIX; 

7EV , TEVR , TEVI : LISTR; 

3RLA5TX , S2LASTY : EXTZRMAL INTEGER; 

TEMPI CM , T EMP I CM: I ME EG E R ; 

G3CET : REAL; 



(* 



XT E REAL 


PROCEDURES 


amc mulct i cm 


S *) 


EXTERNAL 


PROCEDURE 


CRT 1 tl IT; 




EXTREMAL 


PROCEDURE 


CRT ( C : CRTCC ' 


LAND ) ; 


EXTERNAL 


PROCEDURE 


ppmp\M? / v v • T 

JD i VAi ^ \ / i • 1 


NT EGER) ; 


EXTERNAL 


PROCEDURE 


PFOMPTAT { Y : I 


NT EGER ; S 


EXTERNAL 


PROCEDURE 


CLEARSCREEN; 




EXTEF.MAL 


PROCEDURE 


CLF.ARIT ( I : IN 


TECER) ; 


EXTERNAL 


FUNCTION 


GZTCNAR (ONSET 


• r- ■n m t ♦ 

. ^Li Ji. ^ • * 


EXT EE ME L 


PROCEDURE 


GET STRING (VA 


R G : G T R I 



• ^ ti ^ p 



GjMAXLE: 



HER) ; 



1 



rr ) n 

U H 



EXTERNAL FUNCTION YES: BOOLEAN; 

EXTERNAL PROCEDURE WAIT; 

EXTERNAL PROCEDURE NHEAD (S : STRING) ; 

EXTERNAL PROCEDURE INTREAD (VAR K: INTEGER); 

EXTERNAL FUNCTION VALUE (VAR S : STRING ; VAR p : INTEGER) : 

PR. I NG ; MA.XLEN : INTSGI 



EXTERNAL 


PRCC 


EDURE GETRE AL ( VAR 5 


EXTERNAL 


[3] 


PROCEDURE MATRIXIC; 


EXTERNAL 


[20] 


PROCEDURE HELP; 


n A X iA ft Li 


[19] 


PROCEDURE TTRDET; 


^ y rn n p s f y 


[20] 


PROCEDURE TTG-'XB; 


EXTERNAL 


[21] 


PROCEDURE TTRPQR; 


EXTERNAL 


[24] 


PROCEDURE TTRNAA; 


END OF EXTERNAL DECLARATION'S *) 



REAL ; 
R) ; 



P ROC ED U RE P REAXB ; 

(* SAVES 3 WHILE THE INVERSE OF A IS CC ’PUT EC *) 
VAR I , J : INTEGER; 



3 EG IN (* REE, AX 13 * 



IF ICEFLAG THEN 
3 EG IN 

TEMPICM := ICG; 

FOR I : = 1 TC ROND IMS DC 
FOR J : =1 TO ION DO 

TEMPTS [ I ,J] := I C T 3 [ I , J ] ; 

TEMP I ON := ROND IMS; 

END ; 

IC := ICNDIM; 

FOR I : = I TC ICNDIM DC 



FOR 



TC ICNDIM DO 



ICTB [ I , J ] := 0.0; 

FOR I : = I TC ICNDIM DC 
I0T3 [1,1] : = I . C ; 

TTSAX3 ; 

IF IG3FLAG THEN 
BEGIN 

I CM := TEMPICM; 

FOR I : =1 TO RONDIMB DO 
FOR J : = 1 TO ICM DC 

I CT3 [ I , J ] := TEMPTS [ I, J]; 



END; (* PEEAX3 *; 



2 



BEGIN (* MATRIX *) 

CRT I MIT; 

QUITFLAC := FALSE; 
AO I FLAG := FALSE; 
I02FLA3 := FALSE; 
REPEAT 



\ -Jt-ILV f 

: r? ? r\ t • r t *■ r 



GC 



'.ye;d( 

3CTGXY ( C , 3 ) 

WEI TEL:’ ( ' ' 

WEITELM; 

• < n ITELM(' ' 

WEITELM ; 

.’RITELM (' ' 

WRIT ELM; 

WF. I TELM ( ' ' 

WEITELM ; 

'."RITELM ( ' ' 

WEITELM; 

MF I TELE ( ' ' 

WRITELM; 

WRITELM (' ’ 

WRITELM ; 

.T.ITELM; ( ' ' 

WEITELM; 

WRITS ( ' ' : f , ' SELECT OWE 

I L AC I 2 lj t- m i H :j I 
DFG7M 

| I ^ I » p t _ 

IF I 03F LAG THEM 
2 EG I M 



rr i r> 
J , u 



H t ^ 
- f ^ 



r » t? 

- ; L 

r: i ^ 

- f ^ 

3, 'C 



^2*. r r 1 r a ■ i 

PEAL PCLV IOMIAL ROOT GCL"GR ' ; 

F I GEMVECTCFG CF RE.' L 'ATPIX'); 

SOLUTION CF MCMSIMGULAR LI MEAT GQU .7 1 CMC ' ) ; 
DETERM I M AMT CF REAL MATRIX’ ) ; 



r » t 7 — ^ • -p r- ’-i 
iM'/llhCL 



MEAL MATRIX’); 



ED I T I MG AMD DISPLAY CF IMFUT LY 



■ r* — ■ 'T I \ 



D I RLCTICM3 FOR USE <F CALCULATE I ' ) 



, 'QUIT' ) ? 

: ')• 



:f xcwdimb=ioncim them 



/KSST := CKSET+ [ 1 C ' 



M • 

*-**■' L-' / 



HD 



ELSE 



’A'/F' 



] ; 



OKCET := OXCET+ [ ' G 1 , ' Q 1 ] ; 
jlLEc i' z — j wT'Li . . R ( 0 1 \ 3 E T ) } 
^ l \ R C H. Ci E ix } 



BEGIN 

FR0HP7AT ( 1C , 

if yes the:; 

CEITFLAG : 
CLEARCCEEEH ; 

end 



'nn v'mt iz 



'■ E THE EDITED •1DTRK 



Y/" ’ ) 



TFJE; 



ELSE 

CASE 



ELECT CF 
TTRPQR; 
TTF.rir 

THP ' « r . 

TTRDET ; 
PREAXB ; 

HC TEIXIC; 



3 



HELP 



'G' : 

END; 

UNTIL QUITFLAG; 
END. (* MATRIX *) 



(* VERSION 0289 



MODULE CETLI3; 



:rtco:imand = 



ij i. w l! ^ .in — • j *j 1 J i 



(ERASECS , ERASEOL , UP , DOWN , EIGHT , LEFT, LEAPT' 
FCCL r *R , 3CCLCR , RE VI DON , REV I DOF" , IN TENON , I 

p.linxcm, 

3LINRCFF) ; 



'* r 1 

tf: 



T*- 

r 



:cnin 

? r rv v 



1 t ^ * * 
. 7 rp r-t r- r/ c* r m 
; l\ ^ iv L 1 L r i\ ^ J L. i 



PTR = ''INTEGER; 

EP.iCPE RATIC'I = (C0LD3CCT , NARH3C0T , C 0 IIS TAT , 

PUHCUT , RDR III , HOME , S ELDER , E 
DUKREAD , D3 K\ .RITE) ; 

COLOR = (BLACK , BLUE , GREEK , CYAN , RED , HA GELT A , BRCHN , LGR AY 
LGREEH , LCYAM , LRED , L.1AGEKTA , YELLOW /.HI I TE ) ; 
3TRIMG4C = 5TRIHG [40] ; 



i 



Y,L3LU; 



BELL = 07; 

RTM = 



= 



1 ■; . 
J — / 



o ; 



VAR 3BLASTX , GEL AST Y : INTEGER; 

CRT I LIFO : ! RRAY [CRTCCH 1AHD] OF CHAR; 
COLOR I NFG : ARRAY [COLOR] OF INTEGER; 
PREFIXED: ARRAY [CP.TCC MAKE ] OF BOOLEAN; 
EDCUVAL: INTEGER; ( ^ GLC3AL VARIABLE FO 



J. L Ji 



CALLS * } 



ETERNAL FUNCTION GRDCC00 (FUNC : INTEGER ; ?\RH:?TR) 



PROCEDURE 3ICSCALL ( FUNG : CP.HOPERAT ION ; CCH : INTEGER) ; 



VAR RESCRIPT: ARRAY (1.. 5] 
T3ITE: BYTE; 

J: INTEGER; 



IYTE ; 



I IN (* F IOGCALL *) 
F FUNC=CCNOUT THE 



!EGIN 



— r n u , 



HABIT (T3ITE, 8) ; 



DESCRIPT [1] 
DESCRIPT [2] 
DESCRIPT [ 2 ] 
DESCRIPT [4 ] 
RESCRIPT [5] 



= .1 



TBITE; 

I 1 • 
r 

i i . 
r 



EDCSVAL := 3 BROS 8 6 ( 5C , ADDR (DECCRIPT [ 1] ) ) ; 

END ; 

T [7 P'1’ I ^ — (- r- ’ T T « l rntjn"- 

- -- i_ - I — ' — v . l ■. i 1 u uil 

BEGIN 

DECCR I ?T [ 1 ] := 2; 

FOR J : = 2 TO 5 DC 
RESCRIPT [J] := ' 



5 



GO o 



3D0SVAL := 0BDOS8 6 ( 50 , ADDR (DESCRIPT [ 1 ]) ) 
END ; 

END; (* 3 IOSCALL *) 



PROCEDURE CRT I NIT; 

V Tv Ji. O A i • L\ ± 1 *. 1 L « l 

33 GIN (* CRT I NIT *) 



' / 



Cm .170 [LEADIN] : 


= OHR (27); 


CRT I NFC [ERASEOS] 


: = 1 J ' ; 


CRT INFO [ERAS EC L ] 


• — 1 y 1 . 




CRT INFO [RIGHT] := 


i r* > • 

e 


CRTINFC [LEFT] := 


1 D' ; 


CRT INFO [UP] := ’A 


i . 
/ 


CRT INFO [DOWN] : = 


' 3 ' ; 


CRT INFO [TINE] : = 


1 1 • 
f 


CRTINFC [FCCLOR] : 


= "o'; 


CRT INFO [BCG LOR] : 


_ mi. 

“ ^ r 


CRTINFC [ REV I DON] 


: - 1 p 1 } 


CRTINFC [ RE VI DOFF ] 


• — 1 /"«■ 1 • 
• “ 4 / 


CRTINFC [INTENON] 


. __ 1 y 1 . 

* “ L f 


CRTINFC [ I NTENCFF ] 


- _ 1 ■ , t . 

« / 


CRTINFC [’3LI IKON ] 


• — 1 Q 1 • 

« ^ t 


CRT INFO [3LINKOEF] 


• - 14- I . 

• ^ f 


PREFIXED [LEADIN] 


:= FALSE; 


PREFIXED [ERASCOS] 


:= TRUE; 


RRCFIXED [ ERASECL ] 


: = TRU E ; 


PREFIXED [RIGHT] : 


= TRUE; 


PREFIXED [LEFT] := 


’v IE; 


PESEIXED [UP] := T 


RUE; 


PREFIXED [DOWN] := 


1 iA U , 


PREFIXED [TINE] := 


l nLoE ; 


FOR OPCRT: =FCCLOR 


TO DLINKC 


PREFIXED [ CPCRT ] 


:= TRUE; 



END; (* CRT I NIT *) 



PROCEDURE CRT ( C : CRTCC.NMAND ) ; 

DEGIN (* CRT *) 

IE PREFIXED [C] THEN 

3 IOSCALL (CONCUT , ORD (CRT I MEG [LEADIN'] )) ; 
3 IOSCALL (CONOUT , ORC (CRT I MFC [C ] ) ) ; 

END; (* CRT *) 



PROCEDURE GOTOXY (X,Y: INTEGER) ; 
VAR I: INTEGER; 

BEGIN (* GOTCXY *) 

FOR I : =1 TO LOG DC; 

CRT (LEADIN) ; 

FOR I : =1 TO 100 DC; 

2 IOSCALL (CONCUT, CRD ( 1 Y’ ) ) ; 
FOR I : = 1 TO ICC DC; 

C IOSCALL (CONCUT, Y+3 2) ; 

FOR I : =1 TO ICQ DO ; 
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BIOSCALL (CCNOUT, X+22) ; 
SBLASTX := X; 

SBLAETY := Y; 

FOR I : =1 TO 100 DO; 
END; (* GOTOXY *) 



M^n^r I v . T^ r 'r^r i n . C* • r*r-» 
' ■ ' i ^ x . -L ^ ^ — 3 ii, i\ f %j • *j x 



VAR J: INTEGER; 

g. r ^ t ‘] ' * p R C • 1 P T 1 ~ T * ) 

GCTCXY (0 , Y) ; 

JRITE (S) ; 

:rt (erasecl) ; 

FOR J : =1 TO 4 DO 
CRT (TIME) ; 

END; (* PROMPTAT *) 



"G ) ; 



PROCEDURE CLEARGCREEN ; 

•7 . . R U i INTEGER; 

BEGIN (* CLEARGCREEH *) 
GOTOXY (0,0) ; 

CRT (ERAS EC 3) ; 

FOR J : =1 TO 4 DO 
CRT (TIME) ; 

END ; (* CLEARSCREEU *) 



PROCEDURE CLEAEIT ( I : IN r 



i ± U'jU Li 



ry * ’ 



BEGIN (* CLEAR I 
GCTCXY (C , I ) ; 

CRT (ERAS EGG) ; 

FOR J : = 1 TO 4 
CRT (TIME) ; 

END; (* CLEAR IT *) 



DO 



FUNCTION GETCHAR (ONSET 
VAR CM: CHAR; 

OCH: INTEGER; 

GOOD: 300LEAN; 
BEGIN (* GETCHAR *) 
REPEAT 

BIOSCALL (COM I M , 0) ; 
CM := CHR(EDOGVAL) 
OCII : = CRE (CH) ; 



ETOFCHAR) 



SAR: 



* 0 IS A DUMMY VARIABLE 



IF OCH >9 6 THEN 
IF CCH< 122 THEN 
CH := C!IR ( CCH- 2 2 ) 
GOOD := CH IN CK.3ET; 
IF NOT GOOD THEN 
..’RITE (CI1R ( 7 ) ) 

ELSE 



IF CH IN [ ’ 
NRITE (CM) ; 
UNTIL GOOD; 
GETCHAR. := CH; 



r 



mrj'niT 
i. 11 jL% 



*) 



7 



END; (* GETCHAR *) 



PROCEDURE GETSTRIHG ( VAR S: STRING; MAXLEN: 
VAR SI : STRING [1] ; 

STE;iP : STRING; 

CXSET: SET CF CHAR; • 



JIN (* GETSTEINC + ) 



. _ r » ! 



i i i 

J 



] ; 



; 1 



— i i . 

/ 



STEM? : = 1 ' ; 

REPEAT 

IF LENGTH (STENiP) = 0 THEN 
SI [ 1 ] := GETCHAP (OKSET ) 

ELSE 

IF LENGTH (STENP) =MAXLEN THEN 

SI[L] := GETCHAR ( [CHR (RTN) ,CHR ( 
ESP) ] ) 



ELSE 

C 1 



o.fi] 






HAR ( OKSET + [ CHR ( RTN ) , •: 



IF S 1 [ 1 ] I II CKSZT THE II 

CTZHP := CCNCAT (STE.1?,S1) 

iid Li o i_i 

IF 31 [ 1] = CHR (ESP ) THEN 
BEGIN 

CRT (LEFT) ; 

HRITE ( ' ' ) ; 

CRT (LEFT) ; 

DELETE (STEM?, LENGTH (STE.1F) , 1) ; 

UNTIL SI [1] = CNR ( RTN ) ; 

IF LENGTH ( CTEMP ) <> 0 THEN 

S : = C TIE IP 
ELSE 



WRIT 


E(G) ; 


END; (* 


GETSTRIMG 


FUR GTIGN 


YES: 3CC 


BEGIN (* 


YES *) 



END; (* Y 



GETCHAR ( [ ' Y ' , ' M ' ] ) III [ ' Y ' ] ; 

a * \ 



PRCCEDUPE GAIT; 

BEGIN (* WAIT *) 

CLEAR IT ( 5) ; 

PP.CNPTAT ( 10 , ’ PLEASE WA IT... ’ ) ; 
END; (* WAIT *) 



PROCEDURE .JHCAD (A: STRING) 
VAR I: INTEGER; 

BEGIN (* WHEFD *) 

CLEARo „ R — I . ; 

I := (8C-LENGTH (A) ) EIV 



INTEGER) ; 



IHR (SSP) ] ) ; 
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GGTOXY (1,0) ; 

WRITELN (A) ; 

GCTOXY (1,1) ; 

FOR I : = 1 TO LENGTH (A) DC 
IF A [ I ] = * ' THEN 

U'RITE(' ’) 

ELSE 

■:EITE( ) ; 

. , R I T E L l ' ; 

ND; ( * HHEAD* ) 



RCCSDURE i: TREAD (VAR K: INTEGER); 

YPE STRINGS = STRING [6]; 

AR .3: STRINGS; 

PROCEDURE GETISTRIUG (VAR 3: STRINGS ; KAXLEN : INTEGER) ; 
VAR 31: STRING [ 1 ] ; 

STEMF , TTEt'iP : STRINGS; 

ONSET, CKAY3ET: SET OF CHAR; 

’ LAG2 , NFLAG : ECCLEAN; 



* -> - / 



T r • \ ", " ' rn - m rpo r-/ pr. T'M 



PECCDDULS CHEC::K3T; 

VAR I: INTEGER; 

3E3IJI(* CHECKI :IT *) 
TTRlIF := STEV.F; 

DELETE (TVEi-lP, 1,1) ; 

Is! ^ i I bilr) < J i j 

T7 r ' ■ '*! 1 - - TT^r-irr 

. .JT', O 1 • — ^ I— 

P r r.' 



BEGIN 




T1 


= CRD (TTSMP [1] ) 


T2 


= CRD (TTSMP [ 2] ) 


T3 


= CRD (TTEMF [3] ) 


T4 


= CRD (TTENP [4] ) 


15 


= CRD (TTEMP [ 5] ) 



IF T 1 < = 5 1 THEN 
IF 7 1=51 THEN 
IF T 2 < = 5 0 THEM 
IF T2= 50 THEN 
IF T3 < = 5 5 THEM 
IF T3 = 5 5 THEM 
IF T4 < = 54 THEN 
IF 74=54 THEM 
BEGIN 

IF NFLAG THEN 
BEGIN 

IF T 5 < = 5 6 THEN 
FLAC1 := FALSE 
ELSE 

FLAG 1 := TRUE 

END 

ELSE 

BEGIN 



9 



cu o 



IF T5< = 5 ! 

FLAG1 : 
ELSE 

FLAG 1 : 

END 



P'jn 



ELSE 




• 


FLAG1 


: = 


F A L S E 


ELSE 
FLAG i 


. = 


C RUE 


ELS E 
FLAG 1 


- — 


FALSE 


ELSE 

FLAG1 


. _ 


TRUE 


ELSE 

FLAG1 


. _ 


FALSE 


ELSE 

FLAG1 


. _ 


TRUE 


ELSE 
FLJ 31 


. _ 


FALSE 


ELSE 

fl; G1 


. — 


TRUE; 



c:ic ; 

IF FLAG 1= TRUE THEM 
BEGIN 

FLAG I := TRUE ; 

L := LENGTH (STEM?) ; 
FCP. I : = 1 TC L DC 
CRT (LEFT) ; 



I 1 



’.JR I TE ( 1 ’ ) ; 

FOR I : =1 TC L 
CRT (LEFT) ; 



DO 

DO 



WRITE (CUR (3ELL) ) ; 

i?Mn • 

*-*-*L' t 

END; (* CHECK IIIT *) 



£ EG I M ( * GET I STRUG *) 

EAT 

KSST := [ 1 C 1 . . ’ 9 ’ ] ; 

SI := ’ •; 

C |T> T? .( o • - II. 

UFLAG := FALSE; 

REPEAT 

IF LENGTH (STE, IP) = C THEM 
BEGIN 

MAX := MAXLEU-1; 



£ L 2 l — 


FALSE; 


OKAYS ET 


:= 0KSET+ [ 


Sl[l] := 


GETCHAR (0 


IF Cl [ 1 ] 


IN (['+', 


EEGIlJ 





IF SI [!] = ’-' THEN 
MFLAG : = FALSE; 



THEM 
= FALSE 

= TRUE 
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MAX := MAX+1; 

FLAG2 := TRUE; 

END; 

END 

ELSE 

IF LENGTH (STEMP) =MAX THEN 

Sl[i] := GFTCHA? ( [CHE (RT") ,CHR 
(3C?) ] ) 

1 L S I] 

L/ i M 

IF (LENGTH (STE I?) =1) AND FLFG2 
THEN 
SI [I] 

GLEE 
31 [ 1] 

END; 

IF S 1 [ 1 ] IN (CKSET+ ] ) THEN 

(' i HrMn . r* , ' * 7^ n m / r* m ^ f r> ^ R 

e i Ltir . — LU JUii i (i 1 ir / cii 

) 



GETCHAK (CKSET+ [CHR (ESP) ] ) 

GETCHAR (OKSET 4- [CHR ( RTN ) , EHR (EE ? ) 



IF Cl [1] = CHR (DSP) THEN 

pj * T 

O - J 1 u 

r "'7T rn f r ^ C7 r y ' • 

y nr / i » \ . 

,\L\X L I- [ ) f 

CRT (LEFT) ; 

DELETE (STEM?, LENGTH (STEM?) , 1) ; 

^ iJ ; 

UNTIL SI [13 = CHR (RTN); 

CHECK I NT ; 

UNTIL NCT FLAG 1 ; 

S := STE, IF ; 

END; (* GET I STRING *) 

PROCEDURE STRTCINT (VAR 3 : STRINGS ; VAR K: INTEGER) 
CCN3T Z = 48; 

VAR STEM?: STRING'S; 



FLA CP: 

I , L : IHTE' 



' r OLE AN ; 



’RTOINT *) 



GTE IP := 3; 

IF STEMP [!]='-’ THEN 



r r cf V 

jijOlj 

FLAG? := TRUE; 

IF (NCT FLAGP) CR ( STEM? [ 1 3 = ' + ' ) THEN 
DELETE (STEMP, 1,1); 

L := LENGTH (STEMP) ; 

!• := 0 ; 

FOR I : = 1 TC L DO 

K := 10*X+CRD (STEMP [ I ’ ) -Z ; 



IF FLAGP=FALS' 
K * — — K j 
END ; ( * 3TE7C I NT 



’MEM 



31 



3EGIN ( * INTREAD *) 
GET ISTRING(3,6) ; 
STRT0I1IT ( 3 , K) ; 
END; (* INTREPD *) 



: ..THING; VI 



1 . 

- ^ r 



rpt t z*' m ' T "* * r '1 n /’-?•- 

i. L fv-i i v., V . ‘ Li J L ^ 1 4 i\ 

-• s t z i -n 
^ V i « O i 

EMIT = 1MC + 

Z = ; (* OED(C) *) 

v; R A , Y : REAL; 

E , I , J , P2 : INTEGER; 

NEC , NEGEXP , GTE : ECCLEAM ; 
DIGITS: SET OF CHAR; 



-L _ 4 1 _J i. 



FUNCTION FCWP.TEM (EX : INTEGER) 
VAR I: INTEGER; 

T: REAL; 

^ * PC wPTEF * ^ 

: :='o; 

c • i. • Vj j 



REAL ; 



IF ODD (EX) THEN 
CASE I OF 
0 



o 

7: T 
3 : T 
END; 

EX := EX DIV 2; 
I := 1+1; 

C JTIL EX = C ; 
PCNRTEN := T; 

END; (* PCNRTEN *) 



= T * 1 . 0 E 1 ; 

= T * 1 . C E 2 ; 

= T * 1 . C E 4 ; 

= T* i . 0E£ ; 

= 7*1.0216; 

= T*1.0222; 

= jdo4; 

= 7*1.0 E12C; 
= T*1.0C256; 



3CGIN (* VALUE *) 

I := 1; 
r : = C ; 

P 2 := 0; 

GTL := FALSE; 

DIGITS : = [ ’ 0 • . . • 9 ' ] ; 

S := CONCAT (S ( * SAFETY CHARACTER * 

l \ : = 0 ; 

E : — G ; 

NEC := ( 3 [ I ] = ' - ' ) ; 

CHILE S [ I ] = 1 ’ DC 

T : = T-1 • 



IF (S[I]='+*) CR NEC THEN 
I : = I + 1 ; 
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WHILE S [ I ] III DIGITS DO 
BEGIN 

IF S [ I ] = ' 0 ' THEN 
P2 := P2+1 

t-> r r' t? 

BEGIN 

P := P + P 2 + 1 ; • 

22 := 0; 

JTL := TRUE; 

END; 

IF A < L I Is I T THEN 

A := 10*A+CED (5 [ I ] ) -Z 
ELSE 

E := H + l ; 

I := 1+1; 

END; 

IF S [ I ] = ' . ' THEN 



2 : = P + P 2 ; 

I := 1+1; 

IF NOT ( 3 [ I ] 

E> r> r* r ' t 
DJjUIh 

INSERT ( ’ C 
I := 1+1; 



: I GITS) tk: 



»2(I) / 



WHILE E [ I ] = ' 0 ' DO 
BEGIN 

P2 := P2+ 1 ; 

IF :\ < L I * . I T THEN 

•? T? C T * 1 

'J I L . 

A := 1 0 * A + C RD ( S f I ] )-Z; 
E : = L- 1 ; 

END; 

I := I +1 ; 

END ; 

IF GTL THEN 
? := ? + P 2 ; 



NHIL 



[I] IN D I G I 1 



BEGIN 

? := P+1; 

IF A < L I M I T THEN 



BEGIN 

A := 10 * A + CRD ( 5 [ I ] ) - Z ; 

E := E- 1 ; 

END; 

I := I + 1 ; 

END; 

IF 3 [ I ] IN [ ’ E ’ , ’ E’ ] THEN 
BEGIN 

I : = I + 1 ; 

J : = 0 ; 

REGEX? : = (3 [I ] = ’-’) ; 

IF ( S [ I ] = ’ + ’ ) CR NEGEX? THEN 
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I := I +1 ; 

WHILE 3 [I] IN DIGITS DO 
BEGIN 

IF JCLIMIT THEM 

J := 10* J+ORD (S [ I ] ) -Z ; 
I := I +1 ; 

END; 

IE I EG EX? THEM 



ELS E 

C : = E + J ; 

ZMD; 

Y := A; 

IF LEG THEM 
Y : = - Y ; 

IF E<C THEM 

VALUE := Y/PCWRTEM (-E) 
ELSE 

IF EOC THEM 
VALUE := 

Y* PCWRTEM (E) 



VALUE := Y; 

IF ((MOT MEG) AMD (Y=0.0)) THEM 
VALUE := -Y; 

■’* i t r y-i ^ r t i — i i ^ r 

f i' 1 I Li U o l I j — 'J ^ 

T . _ T j_ 1 . 

X • ' / 

S := COPY (3,1, LENGTH ( S ) - I ) ; 

EMD; (* VALUE *) 



PROCEDURE IETREAL (VAR 
VAR SI, 32: STRIMG[1]; 

10] ; 



<7 • CTHTM 



1A Vf • IMr-nr-FT ) . 



Fop* C^PT* t< 
j t i tr • o ; LUiv 



IMTSET,ALPSET,SIGSET,CKSET,MISCSET: SET OF CHAR; 
LOG: INTEGER; 

KAR , CHC ICE , MTCER : INTEGER ; 

FERFLAG, EXPFLAG: EOCLCAM ; 
iZGIM (* GETREAL *) 



ALPS CT 

T'jfpcrn 

SIGSET 



= [ ' A ' . . ' Z ' ] ; 

= [ ' 0 ' . . ' 3 ' ] ; 

= [ ' + ' ,*-•]; 



.1ISCSET := [ ' E ' , ' . 

EXPFLAG := FALSE; 



— II. 



'TEMP : = 

■ i . _ i i . 



PERFLAG := FALSE; 

KAR := 0; 

REPEAT 

LCC := LElIGTH (STEMP) ; 

IF LEKGTH (ETE’.IP) =0 TEEM 
3 EG 1*1 

C K GET : = IISC5ET+ I NTSET + S IGGET ; 
31*1] := GETC’IAR (CKSET-f [EUR (RTM) 

EEC 
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LSE 

IF LE’JGTH (STEMP) =MAXLEtI THEM 
BEGIM 

CKSET := [CHR(RTM)]; 

G 1 [ 1 ] := GSTCHAR (OKSET+ [CHR (BSP) ] ) ; 

pMn 

LGE 

3EGIK 

IF ..OT EM PPL AG THE*.; 

BEGIM 

LOG := LEMGTli (STZ.1P) ; 

G2 [1] := STEMP [LOG] ; 

IF S 2 [ 1 ] IH SIGSET THEM 
CHOICE := 1 
ELSE 

IF 3 2 [ 1 ] III IHTSET THEM 
CHOICE := 2 



ELSE 

T V 



2 2 [ 1 ] = ' . 1 
CHOICE CP 



^ • 

— / 



2 : CKSET 

“3 . v n n m 

• V r> L L i 

:md; 



MISCSET+IMTSET ; 
MISCCET+IMTSET+ [CHR (RTM) ] ; 
M I SCSET + IMTSET+ [CHR ( RTF ) ] ; 



E EG I'M 



r ^ 


• — Lj 


EMGTH (STEMP) ; 


r i 


] : = 


STEMP [LOG] ; 


7 i-i -* 

-L i. 


2 [11 


IM JIGSET THEM 


CHOICE 


: = 1 


ELSE 






IF 


32 [ 


1] IN I MTSET THEN 




CHOI 


CE : = 2 


ELSE 






IF 


S2[ 


1 ] = 1 E 1 THEN 



CHOICE := 
: CHOICE c 



1: OKSET 
2: OKSET 
2: OKSET 



IUT3ET+ [CHR (RTM) ] ; 
S I GSET+ 1 MTS ET ; 



r\m . 

— - ‘ ^ r 



= GETCHAR (CKSET+ [CHR (BSP) ] ) ; 



SMD; 

31 [1] 

L-iL* U f 

F Sl[l] IM (CKSET- [CHR (RTM) ] ) THEM 
BEGIM 

STEMP := COMCAT (STEMP, SI) ; 

IF SI [ 1 1 = ' E ' THEM 



EKPFLAG := TRUE; 



] ; 



LED 

ELSE 
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IF SI [1] = ' . ’ THEM 
BEGIN 

PERFLAG := TRUE; 

MISCSET := MISCSET- ; 
END; 

END 

ELSE 

IF SI [ 1] =CHR (3SF) THEM 
BEGIN 

LCC := LENGTH (STEM,?) ; 

IF STEMP [LCC] = ' . ’ THEN 
PERFLAG := FALSE 
ELSE 

IF STEM? [LOG] = ' E ' THEN 
EXPFLAG := FALSE; 

CRT (LEFT ) ; 

,-JFITE ( ' ’ ) ; 

CRT (LEFT) ; 

DELETE (STEM?, LCC, 1) ; 

^ A I TA • 

*- J *-'« *-*' / 

T 77 ' j S' rri V T r r ^ r* : r rvt 

i L jC j. JjAx r Li-1'.j i > 1 J , Li 

MISCSET := MISC3ET+ [ ' E ' 1 ; 



:f mo 



r 1 ^ 



ERFLAG THEM 



MISCSET := lISCSET+[ 
jMT I L SI [ 1 ] =CHR ( RTM ) ; . 

IF LENGTH (STEM?) <>0 THEN 
ELGIN 

b • — STEilP ; 

NTOER := 23-LENGTH (S) 
IF NTOER> 0 THEN 

.vTITEC NTOER); 

— "D • 

END; (* GETREAL *) 

PROCEDURE SPACEBAR; 

VAR CH: CHAR; 

ELGIN (* SPACEBAR *) 

NR I TELL! ; 



*■ ~ T' T m — i t 

. 1 1 - i i £. ( 



y R ( ' pr.r.M: 



SPACEBAR 



CH := G ETC II A R ( [ 1 
■JRITELN; 

END; (* SPACEBAR *; 



(* THE FOLLOWING ROUTINES ARC FOR COLOR /1CNITCPS 



(* PROCEDURE IN ITCOLCR ; 

VAT SHADE: COLOR; 

I : INTEGER; 

BEGIN 
I • =o • 

E DR SHADE : =3LACH TO ’SHITE DO 
BE DIN 

COLOR I NFC [SHADE j :=I; 
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I : 
END; 
END ; * 



i + i; 



(* 



PECCEDUEE ALTCOLCR (C: COLOR) ; 

VAR J: INTEGER; 

DEC IN 

r' r~ - .- r-v x . 

— - *-X —i v_/ 1\ J f 

"OR J : = 1 TO 100 CO; 

3IC3CALL (CONCUR, CCLCRI NFC [C] ) ; 
FOR J : = 1 TC 100 DO; 



(* 



PROCEDURE PAINT (X, Y /WIDTH , DEPTH 
VAR J: INTEGER; 

3 EG IN 

GCVOXY ( X , Y ) ; 

ALTCOLCR (SHADE) ; 

FOR J : =1 TC DEPTH DO 
BEGIN 

NPITEC ’ ; I DTII ) ; 

GCTCXY (X, Y+J) ; 

Lj Lt 0 f 

FOR J:=l TO 10C DO; 

ERD; *) 



EGER ; SHADE : COLOR 



T - ' \ 

IT, ) 



1C D 



TT r T 

i-i • 
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(* VERSION 0286 *) 

MODULE OVERLAY 1 ; 

(* MODULE LIMEYS *) 

TYPE DCMA I’ll = 1 . .20; 

MATRIX = ARRAY [DCMAIM1 , DCM.AIM1 ] OF REAL; 
LI ST I = ARRAY [COMA I*H1] CF INTEGER; 

LISTS = ARRAY [DCMAIMl] CF REAL; 

EXTERNAL PROCEDURE SPACEBAR; 



PROCEDURE RLUC (MD,M: INTEGER ; VAR KER: INTEGER ; VAR ALU : MATRI X ; V A 

L I ST I ; VAR SCALE : L I STR) ; 

VAR I , IND, I PI, IS , J,K,MN: INTEGER; 

3 IG , EL , PIVOT , ROWNRM , T : REAL ; 

CH: CHAR; 



L\ J * * 



FUNCTION AMAX1 (A, 3: REAL) : REAL; 
BEGIN (* AMAXi *) 

IF A<3 THEN 
AMAXi := 3 

TR r 0 1“* 

4jLiiD 1j 

AH AX] := A; 

END; (* AMAXI *) 



BEGIN •* RLUD *) 

• ' r 

I * ip nr-> ■ y ” 



T 1 = 1 P / T V 



V INFINITY NORM CD RON [I] CF *) 



/ Lj .L - I 

RC.JNRM := 0.0; 
FOR J : =1 TO N N D 



- i 

T? . 



(* LU DEC 
(* EXPLIC 
IS : = i 





- 


A, 


MAXI (RC.JN 


‘OV/NRM = 


0 


. 0 


THEM 




' i, T MRM : 


= 


1 


r\ . 
• o , 




•E[I] = 


= 


-L. 


. O/RC, 


;l;rm 


JMPCGIT 


I 


ON 


3Y GJ 


■n r t r* r» 
1UW J 


T ROW 


T 

X 


» t m 
- . i 


ercha: 


ire i 
i o LL> 


. TO MM 




DO 







WRITE 
IND := I; 

IF IONN THEN 
BEGIN 

DIG : = C . 0 ; 

FOR K : = I TO NN DC 
BEGIN 

T := SCALD [ K ] *ADS (ALU [X, I ] ) ; 
IF T> 3 1 G THEN 
DEC IN 

IND := K; 
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c n v 



»~3 Crl 



BIG : = T ; 

END ; 

END ; 

IF 3IGOO.O THEN 
BEGIN 

IF INDOI THEN 
BEGIN 

FOB J : = I TO NN DC 
EEC IN 

T := ALU [ IND , J ] ; 

ALU [IND, J] := ALU [I , J] ; 

ALU [ I , J ] := T; 

el;d ; 

SCALE [IND] := SCALE [ I ] ; 

IS : = - IS; 

END; 

I PI := I +1 ; 

PIVOT := ALU [I, I] ; 

FOE I< : = I P 1 TC NN DC 
BEGIN 

EL : = - ALU [ K , I ] /PIVOT ; 

ALU [ K , I ] := -EL; 

IF ELOC.C TEEN 

FOR J : = I PI TC NN DC 

ALU [ K , J ] := ALU [K, J] +EL*ALU [.I , J] ; 

rp * 1 H • 



TP • 

i_ i_> f 

END; 

:? ALU f I , I ] = C . 0 

r . = a . 

r * f r t 1 •— T ■ T o • 

• •* L ± i * - 1 * J f 



END; 

JN [NN] := IS; 
KER := C; 

END; (* RLUD *) 



PROCEDURE EFDS ( ND , N : INTEGER; VA R KER: INTEGER; VAR ALU: NAT?.: 

L I ST I ; VAR X : LISTR) ; 



X ; V R J 1 



VAR I , V , KP1 , L,L?1 ,N.U ,NN: I N'TEGEE ; 

Z: REAL; 

CL: CHAR; 

’BEGIN (* REES *) 

NN := N ; 

II [NN ] =0 THEN 
GIN 

Y n • = 7 . 

r j . . . - / 

NRITELN; 

’..•P.ITELN ( ' IN RF3S , THE TRIANGULAR FACTOR U OF A ' ); 



in t n 



TELN ( ' IS SINGULAR. A UNICUE SOLUTION DOES ') 

'OTrprtrM/ 1 »T^m ” v T I \ . 

>- 1 1 Lull \ I .J i\ 1 iD* 1 • ) f 

> L “i ' ^ i_. u i * 1 «. / 

i-i A i 7 J 

TT ' 7 ^ • 
u.Jl; f 

V [DO • - n . 

— * • “ / 
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NM1 := NN-1; 

IF NM lOO THEM 
BEGIN 

(* SOLVE LY =3 ( FOR HARD SUBSTITUTION) *) 

FOR L : =1 TO NM1 DC 
3 EC I :j 

,;eite ' 

v- . _ t v t r r l . 

A . - o l. !_ UJ J , 

Z : = X [ K ’ ; 

X[K] := X [ L ] ; 

- • i L ] • = / 

LF1 := L+l ; 

FOE K : =LF 1 TO MM DC 

X [ K ] := X [K] -ALU [K, L] *Z ; 

END; 

(* SOLVE UX= Y (BACKWARD SUBSTITUTION) *) 
FOR I : =1 TO MM 1 DO 
BEGIN 

' T - 1 TE ('...'); 



V * i\ 

K := MM- I ; 
— K 



X i ? 1 : = K + 1 ; 

X[NP1] := X[XP1]/ALU[KP1,KPI] ; 
Z : = -X [KPI] ; 



X [ L ] := X [ L ] +ALU [ L , XF 1 ] * Z ; 



X[L] := X [11/ALU [1,1] ; 



END; ( 



* ?. r r * 



> x ) 



k /“' P' T? * A 

U 7L. J • 
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(* VERSION 0235 *) 

MODULE OVERLAYS; 

CONST MAXDE3 = 20; 

MAXDEGP 1=21; 

TYPE DOMAIN! = 1 . . 2f ; 

CO‘!R?l = l. . 'AXDEGFI; 
DIETER A=A?.EAY [DOMI.Rl] 



IF REAL; 



MATE I X = ARRAY [DOMAIN! , DOMAIN!] CF REAL; 



'I = ARRAY [DOMAIN!] 






LISTR = ARRAY [DOMAIN!] CF REAL; 

CRTCCM 1AND = (ERASE CS , ERAS EC L , UR , DC. /N , RIGHT , 
FCCLOR , DCCLOR , REV I DOM , REVI DOFI 
3L I MKOH , 

3LINK0FF) ; 

SETCFCHAR = SET OF CHAR; 



, 1ST 



-p 7" M 



1 U 



1 n* ~ n - 

J ;j £\ / 



CP MOP ERR riCN = (CCLD3GCT, JAR 13CC T , CONST AT , CO 

RUNOUT , RDRIM , HOME , 3ELDSX , SET 
D5KEEAD , DSR /RITE) ; 

3TEING4C = STR I MG [40]; 

STRINGS = STRING [A] ; 

STR I MGS C = STRING [SO]; 



MIM , CONCUR , L: 

^ i-- * 



•\R AC I FLAG , I CDFLAG : EXTERNAL TOOLE?' ; 

IOM, I CMD I M , ROMD I MB : EXTERNAL INTEGER; 

ICT' , TVEC : EXTERNAL MATRIX; 

1 373, XX: EXTERNAL ’MATRIX; 

TEV , 73VR , TEVI : EXTERNAL LISTS; 



ZZDET: E 


XTERNAL RE 


A L ; 


xternal 


PRCGEDURES 


AND FUNCTIONS *) 


external 


PRCGEDURE 


CFTINIT; 


EXTERNAL 


PROCEDURE 


CRT(C:CRTCGM ’.AND) ; 


EXTERNAL 


PROCEDURE 


GCTOXY ( X , Y : INTEGER) ; 


EXTERNAL 


PROCEDURE 


PT.OMPTAT (Y: 1'ITEGEF ;S:STRING) ; 


EXTERNAL 


PROCEDURE 


CLEARSCREEN : 


EXTERNAL 


PROCEDURE 


CLEARIT ( I : INTFGER) ; 


EXTERNAL 


FUNCTION 


GETCMAR (CKSET : SETOFCHAR) : CHAR; 


EXTERNAL 


PROCEDURE 


GETSVRING (’’AR S : STD I NC ; MAX LEM : 


external 


FUNCTION 


YES: 30CLEAN; 


EXTERNAL 


PROCEDURE 


".M- IT; 



t-jtege:;) ; 



21 



O L‘3 



EXTERNAL PROCEDURE '.JHEAD ; 

EXTERNAL PROCEDURE INTREAD (VAR K: INTEGER); 



EXTERNAL 


FUNCTION 


VALUE (VA 


P C • C m PT M ^ • 

J • ED x I\ X y O f 


VAR 


P: INTEGER): REAL; 


EXTERNAL 


PROCEDURE 


GETREAL 


(U*\D C • 4 T ^ T ' 


? ^ • * ! 7. 
t * 


^ ' r TTTT. T ' T ’T’ TP ^ T7 \ • 

. . U i-i l. . i i! * vl u i f 


EXTERNAL 


PFOCEDURE 


SPACED A 


R; 







PROCEDURE G2TCCLUMN (NRCw, J : INTEGER; VAR A : MATRIX) ; 
VAR 1 1 , K , L L , F 4 : INTEGER; 

SREEL: STRING; 

FLAG: BOOLEAN; 

PROCEDURE COLD I SPLAY ; 

VAR IC: INTEGER; 

BEGIN (* COLDI SPLAY *) 

CLEAT IT ( 2) ; 

GCTOXY (C , 2) ; 



BEGIN 

X := IC+1; 
GCTOXY (0 , K ) ; 



■ T r> z’ M ' r i t f' < i t i i - i \ . 
-•M c\ t , , , J , J- ) , 



GCTOXY (9, X) ; 
NF.ITE (A [IC, J] ) ; 






PROMPT AT (22, 1 IS THIS CORRECT? i/ 
END; (* COLD 1 3FLAY *) 

BEGIN (*. GETCOLUMN *) 

FLAG : = iLS L ; 

REPEAT 

CLEAR IT (2); 

COLD I SPLAY; 

IF NOT YES TNG” 

BEGIN 

PROMPT AT (22 , ’ IF NO CHANGI 
FOR LL : =1 TO NRCN EC 
BEGIN 

K := LL+1 ; 

GCTOXY (9 , K) ; 

IF A [ LL , J ] >=0 THEN 
UNITE (' '); 



r> v/v ' i • 

• / -y it 



! FE33 RETURN ' ) ; 



CUD 

:lge 

F L A. E 



GETREAL (3FEEL, 21) ; 

IP LEN'GTH (CREEL) >0 THEN 

A[LL,J] := VALUE (SREEL, P4) ; 

IND; 

:= TRUE; 
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UNTIL FLAG; 

END; (* GETCOLUMN *) 

PROCEDURE GETPOLGGF (GROW: INTEGER; VAR L : LISTRP1 ) ; 
VAR I , J , NEW I , NRCD IV2 , P4 : INTEGER; 

SREEL: STRING; 

FLAG: BOOLEAN; 



PROCEDURE CC LD I G PLAY ; 

VAR I: INTEGER ; 

BEGIN (* CC LD I SPLAY *) 

CLEARIT (2) ; 

GCTOXY (0,2) ; 

NRCDIV2 := NRC'W DIV 2; 

FOR I : = 1 TO NRCDIV2 DO 
BEGIN 

GOTGXY (0,2*1) ; 

WRITE ( ' COEF CF X'); 

CRT (UP) ; 

JV ITE (NRC.7- I ) : 

CRT (DONE) ; 

TRITE ('='); 

CCTCXY (12,2*1) ; 

WRITE (L [ I ] ) ; 

END; 

FOR I : =NR r 'DIV2- L l TO NR AW DC 
BEGIN 

JEWI := I-NRCDIV2; 

GCTOXY (4 ) , 2*NEWI ) ; 

WRITE ( ' CC CF CF X' ) ; 

CRT (UP) ; 

WRITE (N ROLL- 1) ; 

CRT (DOWN) ; 

WRITE (’ = ' ) ; 

GCTOXY (52,2* NE'JI ) ; 

WRITE (L ( I ] ) ; 

END; 

PI X1PTAT (23 , ' IE THIS CORRECT? 
END; (* COLDIGPLAY *) 



Y/N ') 



BEGIN (* GETPCLCCF *) 

FLAG := FALSE; 

REPEAT 

CL SAP. IT (2) ; 

COLD I SPLAY; 

IF NOT YES TNEN 
BEGIN 

PRCNPTAT (22 , ' IF NO CHANGE - FREES RETURN 
FOR I : =1 TO NRODIV2 DC 
BEGIN 



ICR 



:y (12, 



) -k r 



) 



IF L [ I ] >=C THEN 
WRITE ( ' ’); 

CREEL := ’ 

GET REAL (CREEL, 21) ; 
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IF LENGTH (SREEL) <>C THEN 
L [ I ] := VALUE ( SREEL , P4 ) ; 

END; 

FOR I : =NRODI V2 + 1 TC 'IROW DO 
BEGIN 

MEN I := I-NRCDIV2; 

GOTCXY ( 52 , 2 **N E . ; I ) ; 

IF L [ I ] > = 0 THEN 



GETREAL (SREEL, 21) ; 

IF LENGTH (SREEL) <>0 THEN 
L [ I ] : = VALUE (SREEL, P4) ; 

END; 

END 



U 



HID 



ELSE 

FLAG := TRUE; 
TIL FLAG; 

(* SET P CLCCF *) 



prCCEDLT.S PMEN (1:1 t-JTEGER ; C ; CHAR ; £ : STRING) ; 
BEGIN (* PMEN *) 

CTCXY (0,1) ; 

. TITELN (C : I , ' ' : 2 , S) ; 

END; (* FI'S’! *) 

PROCEDURE HEAD (A: STRING; J : INTEGER) ; 

VAR I: INTEGER; 

BEGIN (* HEAD *) 

I := (CO-LENGTH (A.) ) DIV 2; 

GCTOXY ( I , J) ; 

NRITELN (A) ; 

END; (* HEAD *) 

PROCEDURE MATRIXIC; 

UA pi OH T N’-’r A ^ • *D ''D r ' r ' . 

V .U\ U i i - L) Hv j . • Li Ej : . . f 

CHOICE: CAR ; 

ONSET : SETCFCHAR; 



PROCEDURE EDIT; 

V7.il EDOUIT: 3CCLEAN; 
EDCHO ICE : CHAR; 



PROCEDURE RENTER (A0F3 :CHAR) ; 

VAR I , J , COLD IN, ROND IM: INTEGER; 

NONANE: SPRING [3] ; 

DEGIN (* AENTER *) 

CLEAR IT ( 1) ; 

CASE ACRE OF 



AC I FLAG := TRUE; 
REPEAT 
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PROMPTAT ( 2 ,’ ENTER DIMENSION OF SQUARE MATRIX? ' 
INTREAD ( ICNDIM) ; 

UNTIL ( ( IONDIM>Q) AND ( ICND IM< 2 1 ) ) ; 

RCWDIM := ICNDIM.; 

COLDIM := ICNDIM ; 

MAMA ME : = 'A' ; 

•’fp • 

BEGIN 

1 02 FLAG := TRUE; 

REPEAT 

PRC 1PTAT ( 2 , ' ENTER RC 7 DIMENSION OF G ? '); 

INTREAD (ROND IMS) ; 

UNTIL ( (RONOIMOO) AND ( ROND IM3< 2 ’ ) ) ; 

REPEAT 

CLEAFIT ( I) ; 



PROMPTAT ( 2 ENTER THE COLUMN DIM 
INTREAD ( I CM ) ; 

UNTIL ( (10 M > 0 ) AND (IOI-K23)); 
COLDIM : = 10 1; 

ROND I M := F.CWCIME; 

M A NAME : = ' 2 ' ; 



n •? i 



EMI 



FOR I : =i TO RCWDIM DC 
FOR J : =1 TO COLDIM DO 
CASE ACRE OF 

'A': IOTA [ I , J ] := 

' 3 ' : I CT3 [ I , J ] : = 0 . C ; 

T 7 »in . 



0 



JL jL. 

CLEARIT ( 1) ; 

NRITZLM (' ENTER COLUMN ' , J , ' OF MATRI] 
CASE ACR3 OF 

GETCCLUMN (ROWDIM, J, IOTA) ; 
CETCOLUMN (RCWDIM, J, ICT3) ; 



MAN AM E) ; 



'A' 

i -} i 



END; 

END; 

CLEARIT ( I) ; 
END; (* AENTER * 



PROCEDURE AED I T ( ACR3 : CHAR) ; 

VAR AC HO ICE: CHAR; 

AQU I T , EFLAG : BOOLEAN; 

I , J, COLDIM, RCWDIM: INTEGER; 

1ANAME: STRING [3]; 



! EG I N (* AEDIT *) 



r* r ^ t-i 

1 A * 



'O n ' 

jl 

: BEGIM 

EFLAG : 
MANAGE 



a p T r ? r 7i ^ • 

..O 1 L uMa; j 

-■ ' A 1 ; 

T O'lD T M • 
i J li i f 
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n r; 



END 

IN 
li L \ 




COLDIM : = IONDIM; 
END; 

BEGIN 

EFLAG := IC3FLAG; 
NAMAME 
RCNDIN 
COLD IN 
END; 



= RCvIDINB ; 
- ION; 



F NOT L> AG THEN 
\ENTER ( ACR3 ) 

LSE 

-:egin 

ACUIT := FALSE; 

REPEAT 

CLEAR IT ( 1 ) ; 

PMEH (2, 1 V , ’ EDIT CURRENT iATRIX ' ) ; 

A NEN I1ATKIX N) ; 



» r>. r ] t m » \ • 
~ Jr 



PUEN (4 , ’ E ’ 

5' : ’C 

GCTCXY (2,8) ; 

MRITL ( ' SELECT CUE : ’ ; ; 

CKSET := [ \A ' , ' 13 1 , ' Q ' ] ; 

ACHCICE := GETCHAR (OKS ET ) ; 

CASE ACHOICE OF 

pro j • — *" 'tj p ^ Q f r“ 

b ecu: 

CLEAR IT (1) ; 

rEITELN (' ENTER COLUMN ',1,' CF !. 
CASE ACRE CF 

' A ' : OETCCLUUN { RC.7D IH , I , ICTA ) ; 
’B’: GETCOLUMN ( RC'.'D YA ,1 , I GTE ) ; 
EMD ; 

E N D ; 

'V: BEG IE 

CLEARIT ( 1) ; 

GCTOXY ( C , 2) ; 

trite ('dc ecu tamt to erase cufre:.” 

'■? Y/II ' ) ; 

CRT (ERA 3ECL ) ; 



r* T V ! 



i A X A 



:RT (TIME] 



IF YES T 



t ! ’T 1 V 



i r ' 



A ENTER (ACR3) ; 

EMD; 

AC (JIT := TRUE; 



EMC ; 

UNTIL AQUIT; 

:md; 

(* AEDIT *) 



I.: I -A 



(* EDIT *) 

c ' i~- l. j_. -j ; 

D( 'MATRIX EDITOR' , C); 
U IT := FALSE; 
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REPEAT 

CLEARIT (I); 

PflEN ( 2 , ' A 1 , ' EDIT MATRIX A'); 
?MEN (4, '3' , ’EDIT MATRIX B’); 
PMEtJ ( 6 , ’ Q’ , ’QUIT' ) ; 

CGTOXY (2,3) ; 

•,’RITE ( ' SELECT CEE •: ' ) ; 

O . _ r I -v I 1 -> | I IT. 

^ . vO i-j i . - L t ^ J r 

*> w 4 1 Vv ± « ~ jj ill i. ^ i"l c ( C i\ >. 

CAGE ECCIIC ICE Cr 
'A' : ADD IT ( 'A' ) ; 

’3' : AED IT ( ' B ' ) ; 

’ C : EDQUIT := TRUE; 

END ; 

UNTIL ECQUIT; 

EIID; (* EDIT *) 



’) ; 



PROCEDURE DISPLAY: 

VAP BISCUIT: 3C0LEAN; 
OISCHOICS: CHAR; 



PROCEDURE DDICFLA (ACRE : CHAR) ; 
v; 1 I , J , COLD IM , ECvJD IM : INTEGER; 
CH: CHAD; 

CD LAG : 3 OC L C A N ; 

‘‘.AM A; IE : STRING [3] ; 

SEC IP (* DDIGPLA *) 

CLEARIT ( 1 ) ; 

CASE ACRE CE 
'A': BECI’I 

CFLAG : 

MANAME 
COLD III 
ROWS I M 
END; 

' B ' : 3 EG IN 

CFLAG : 

mama: ie 

COLD IM 
ROND K1 
END ; 



^ T — - - ^ . 

i V 1 j. f 

1 -\ ' . 

• - / 

IOMDIM; 
I CUD III ; 



= I GET LAG; 

_ i p t . 

= ICC; 

= RCwDIME; 



ID NOT CFLAG THEN 
BEGIN 



Nr, r 



'0 X Y ( 0 , 3 



TRITE ( 1 THEE r 
CRT (ERAS EC L) 
FDR I : = 1 TO 
SPACEBAR; 

MD 



IS NC MATRIX ’ , MAMAM! 
.00 DO CRT (TIME); 



el: 



3EGII: 

NRITELN ( ' THE COLUMNS OF MATRIX ' , MA 
FOR J: =1 TO COLDI 1 DC 
SCSI N 
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WRITELN (' COLUMN ',J,' OF ' , MANAHE , ' IS 

FOR I : =1 TO ROWDIM DO 

BEGIN 

WRITE (HA NAME , ' I J ; 

CASE A0R3 OF 

' A ' : NR ITELM ( I OTA [ I , J ] ) ; 

WRITELN ( ICTE ( I , J] ) ; 



END; ( 



• > • 

t-i | t r> . 

f 

END ; 

S ? ACES AR ; 
END; 

END; 

* DDISPLA *) 



PROCEDURE HCCPY ( AORB : CHAR ) ; 
VAR CH: CHAR; 

CFLAG: BOOLEAN; 

N ANA: IE: STRING [3]; 

CNUH , I : 

BEGIN (* HCCPY *) 

•CASi 



..4 1 Jiji\ , 



)U r.UhC Ui: 

BEGIN 



CFLAG := AO I FLAG; 
MANAME := 'A'; 
CNUH := 3; 

END ; 

12 Hi C l A 

CFLAG := I OBELI l; 
1ANAHE := ’3’; 

CMUN := 4; 

END ; 



CLEAR IT ( 1) ; 

IF NOT CFLAG THEN 
BEGIN 

GCTOXY (3,3) ; 
WRITE (’THERE 
« 3 T ( iFASEol) } 

OR I : 

S P AC EC AR ; 



'JC AV 



’ , MA 



FOR I : = 1 TO ICG DO CRT(TIIE); 



M T » V 37 rr » 



) ; 



EL3E 

HARDCOPY (CHUM) ; 



MUD; (* HOC 



) 



3EGIU C 



» I S PI 



) 



DISQUIT := FALSE; 
REPEAT 

CLEAR SCREEN ; 



HEAD ( ’ M7 


TRI 


:: display’ 


C', \ . 

r ^ ) f 




PNEN ( 2 , ’ 


A ’ , 


’ SCREEN 01 


SPLAY OF A’) 


r 


P ’BN ( 4 , ’ 


O , 


’SCREEN D I 


SPLAY OF 3’) 


r 


PH EM (6 , ’ 


» 

^ t 


’ HARDCOPY 


OUTPUT OF A’ 


) ; 


P.iEN (S , ’ 


E ’ , 


’ HARDCOPY 


OUTPUT OF 2 ’ 


) ; 


P 1BN ( IC , 


' Q ’ 


, ’QUIT’ ) ; 
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GOTGXY (2, 12) ; 

WRITE ( ' SELECT ONE : ' ) ; 

CKSET : = [ ' A ' . . ' D ' , ' Q ' ] /’ 
DTSCHCICE := GZTCHAE ( OKSST ) ; 
CASE DISCHOICE CF 

DDISPLA { ' A ' ) ; 

DEI .'FLA ( ’ N ’•) ; 

-• ''QSV ; 

HCCFY ( ' N ' ) ; 

FISCUIT := TRUE; 

END; 

UNTIL BISCUIT; 

END; (* DISPLAY *) 



2EGIN (* /1ATEIXIO *) 
CLEARSCREEN ; 

QU ITFLAC := FALSE; 
REFEAT 

CLEARSCREEN; 



HEAD ( ’ 


i;;?UT 


\7 T R 


IX 


- ' r- i pi \ 


PNEN ( 2 


1-1 1 

t i * / 


77 T ,T1 

iLUl i 


i; 


TRICES ' 


PMEN ( 4 


1 -7 1 1 
/ t 


^ ^ 


Eli 


OR HARD 


P 'lEN ( 6 


1 r\ 1 1 

t V / 


CUIT 


' ) ; 




GOTO XY 


(2,8) ; 








WRITE ( 


'SELECT ON 


T~» . 

u 


' ) ; 


ONSET 


:= ['A 


i i n 

/ ° 


f 1 

/ 


n ' 1 • 

*. J / 


CHOICE 


: = CE 


tct; 


R (C 


XSET) ; 


CASE C 


HCICC 


TP 
KJ c 






’ A ’ : 


EDIT; 








1 T 1 

LJ 


DISPL 


AY; 






’C : 


QUITFLAG 


: - 


TRUE; 


• 

— * - < — ' t 










TIL Q'JITFLAG 








(* MATRIX 10 


*) 







PROCEDURE HARDCOPY (HNUN: INTEGER) ; 
VAR F: TEXT; 

• CH: CHAR; 

PRELAC- : BOOLEAN; 

RESULT: INTEGER; 

PROCEDURE SET PR I NT; 

VAR CH: CHAR; 

FTRIES: INTEGER; 



BEGIN (* SETPRINT *) 

PRFLAG := FALSE; 

FTRIES := 0; 

REPEAT 

ASSIGN (F, ’ LST: ’ ) ; 

RE. .’RITE (F) ; 

IF IORESULT=255 THEN 
BEGIN 

FTRIES := FTRIES+1; 
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IF FTRIES<=2 THEN 
BEGIN 

NRITELN ( ' PUT PRINTER ON LINE '); 
SPACEBAR; 

END; 

END 



< ^ c. 

? R F L . ' 



RUE; 



UNTIL PR FLAG DR 
END; (* GETFP I NT ’ 



TR I E3> 2 ) ; 



PROCEDURE PRSAXB ; 

V.-M\ i . f J ! I ; T E ‘o E R ; 

3E .1 1 N (* PRSAXB *) 

FOR J : = 1 TC IOM DC 
BEGIN 

NRITELN (F) ; 

■7RITELN (F, 'THE SOLUTION FOR COLIN IN ' , J, • IS '); 
FOR I : = 1 TC IONDIM DC 
u . j G i N 

NRITELN ( F , ' [ ' , XX [ I , J ] , ' ] ' ) ; 

END; 

END; 



* 



PS AXE *) 



PE DC 



ECU RE PR 


MATRIX ( 


D I " , I , J : 


INTXGE 


MANAME : 


C* m r- -T X - r* r 
o i .JlJu [ 


N ( * PE '’A 


TRIX *) 


GE ACE.3 


OF 



BEGIN 
LANAME 
DIM := IONDIM; 
EN’D ; 

BEGIN 

MANAME := ' B'; 

DIM 



= ION; 



Cii 



NO 



ij _>( D ; 

NRITELN (F) ; 

NRITELN (F, 'THE COLUMNS CF 
NRITELN (F) ; 

J : = i TO DIM DC 



:IX ', MANAME,' ARE: '); 



t- L\ 



.J O 1 1.1 

NRITELN (F) ; 

,'RITELN (F, ' COLUMN ',J,' CF ', MANAME,' IS :’); 
CASE A0R3 OF 

'A': FOR I : = 1 TO IONDIM DO 

NRITELN (F, ' A [ ' , I , ' , ' , J , ' ] = ' , IOTA[I,J] ) ; 
'3': FOR I : = 1 TC EC 'C I IE DO 

NRITELN (F, ' E [ ' , I , ' , ' , J , ' ] = ' , I C ? 3 [ I , J ] ) ; 

END ; 



END; (* PRMATRIX *) 
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^ r* ' » D r ^ V ^rr*TT7^ f Mnr T I N . 

/ . i. I_J ». • t\ J \J*jH (, u H IJ _/»J i kJ j f 



PROCEDURE PRRNAA ; 

VAR I, J: INTEGER; 

3 EG I N (* PPRNAA *) 

TRITELM (F) ; 

TRITELM (F , ’ THE ARRAY - 
.IF. I TEEN ( F ) ; 

FOR I : =1 TO T CNCI I DC 

TRITELM ( F , • [ • , TEVR [ I ] , ’ , ’ ,TEVI [I] , ’ ] ’ ) ; 
TRITELM (F) ; 

TRITELM (F, ' THE COMPLEX EIGENVECTORS ARE '); 
.■JR I TEL N (F) ; 

J : = 1 ; 

REPEAT 

IF 'TEVI [ J ] < >0 . 0 THEN 
BEGIN 

TRITELM ( F ) ; 

VRITELW (F, 'VECTOR ' , J , ' IS '); 

FOR I : = 1 TO ICMDIM DC 



( F , 



I \ l 



I t ! f I -*■ t -J r t -riT/T^r T 7 1 \ - 

• / v t L /-/ J — t ± J l- 1 "-' [ ~ t \ j r 



TR I TELM (?,',', TV EC [ I , J + 1 ] , ' ] ' ) ; 

FDD; 

TRITELM (F) ; 

tritelm (f, 'vector ' , j+i, ' has components ’); 



C *- N ^ 1 1 L 

<7 r 1 t - t 

i -1 X - J 



» ? - r * 



TRITE ( F , ' VECTG R ' , J + l , ' [ ' , 

TRITELM (F, ' , ' , - TYEC [ I , J + l] ,']'); 



, TV EC [ I , J ] j ; 



• — wv r <- f 

T 7 > T» t rr 

ELSE 

BEGIN 

TRITELM (F) ; 

TRITELM (F VECTOR ’ , J , * HIS COMPONENTS ' ) ; 

FOR I : =1 TO IONDIM DC 
SEC IN 

TRITE (F, 'VECTOR' ,J,'[', !,'] = [' , TVEC [ I , J ] ; 

TRITELM (F , ',',0.0,']'); 

END; 

J := J+l; 

UNTIL J> I GUDIM; 

END; (* PPRLIAA *) 



PROCEDURE PRRDET ; 

BEGIN (* PRRDET *) 

:RITELN (F) ; 

NRITCLN(F, 'THE DETERMINANT OF MATRIX A IS 1 , GSDET) ; 
END; (* PRRDERT *) 



CIN (* HARDCOPY *) 
GET PR I NT; 

IF PRFLAG THEN 
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BEG I N 

CASE HNUM OF 
1: PRRNAA; 

2: PRSAXB ; 

) ; 
) ; 



END; 

END; (* HARDCOPY *) 



BCD END . 



2: PRMATR IX ( ' A ' 
4: PRMATRTX ( ' B ' 
5: PF.RDET; 
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(* VERSION 0235 *) 



MODULE OVERLAYS; 

(* EIGEIJHCR 1 CDULE *) 



i: 


!i = " 


O • 












1 •* 


r — \ "E C T 

A. — K K i\ 


Y [DC !A 


I 


Nl, 


D G l 


1 A 


INI] OF 


T I 


= A FRAY 


[DOM a: 


N 


1 ] 


G F 


T 

J. 


NTEGER ; 


TR 


= ARRAY 


CJ 

O 

? 

KH 


- 


1 ] 


C j 


R 


EAL; 


'HR' 


] A L PRCC 


EDURE 


o 


PAG 


E 37 


oR 


f 



SAL; 



PROCEDURE :>CR2 (MM, N, LOO/, IGH : INTEGER ; VAR H : MATRIX ; VAR NF,NI 
VAR 2 : MATRIX ; VAR I ERR : INTEGER) ; 

CONST , iACHEP = 1.0E-14; 

V.-iR K , E 'M2 , I , 1 1 , I TS , J , T , L , i i , i P 2 , M A : INTEGER; 

NORM , P , Q , R , RA , S , SA , T , TEMP 1 , TGI iP2 , VR , VI , W , X , Y , Z R , Z I , 2 Z : 
F1AG1 , FLAG 2 , FLAG 3 , FLAG 4 , NOT LAS , TESTER IT : BOOLEAN ; 

CM: CHAR; 

FUNCTION SIGN (E,F: REAL) : REAL; 

BEGIN (* SIGN *) 

T T G' / r> n'jr-iT 

a \ c. x 1 1 — i'j 

SIGN := - A 3 3 ( E ) 

E L E 

SIGN := ABS(E); 

END; (* SIGN *) 

FUNCTION NINO (II, 12: INTEGER) : INTEGER; 

BEGIN (* MI NO *) 

T F Il< 1 2 THEM 
MI NO := 11 
ELSE 

MI NO := 12; 

END; (* MI NO *) 



LT 



RLN.L; 



PROCEDURE CCMDIV (XR , XI , YE , YI : F SAL ;VA 
VAR D , H : REAL; 

BEGIN (* 20MDIV *) 

IF A3S ( YR) <AES (YI ) THEN 
BEGIN 

H := YR/YI; 

D := YI+H*YR; 

ZR := (XR*H+XI)/D; 

21 := ( X I * H - X R ) / D ; 

END 
E L A E 
BEGIN 

n := YI/YR ; 

D : = YR+H* Y I ; 



2 R, 21 : REAL) ; 
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ZR := ( XR + H* X I ) /D ; 

ZI := (XI-H*XR)/D; 
END ; 

END; (* COMDIV *) 



PROCEDURE I QRCt-Tl ; 
VAR I: INTEGER; 



PRCC 


L j 


DU 


R 


E 


L 


COK 


VAR 


r 




I 


* T -'I 
14 i 




CER 


3BGI 


L4 


(* 




LC 


c 


K * 


.1 


: 


= 


E 


ZJM 


2 


+ 1; 


FL 


A 


G 2 




: = 




FAL 



REPEAT 

M := 11-1; 

ZZ := H [M , M] ; 



7 7 ” • 

L — Zi l> , 

(R*c-.j)/n[-i+l / 'i]+H[:: / :i+l] ; 
H[ 1+1, I+1]-ZZ-R-S; 

?! [M+2 , U+ 1 ] ; 

ASS (?) +A3S (C) + ADS (R) ; 

P/3; 

C/2 ; 

E/S ; 

I? ;>:=L t:ien 
FLAG 2 := TRUE 

ELSE 



TCE-I := ABE ( H [ M , U - 1 ] ) * ( AE 3 ( C ) 
TEMP 2 := ABS (H [M-l ,M-i] ) +AE3 (Z 
TE I P 2 := 1ACHEP* ABS ( F ) * TEMP 2 ; 
IF TEMPI <= TEMP2 THEN 
FLAG2 := TRUE; 

END; 

wEITS 

.J i.\ i Ii-j l b(i'j2 = i tGi- J 

IP 2 := .1+2; 

FOR I : = IF 2 TC EU DC 
BEGIN 

*1 [ I / I - 2 ] := 3.C; 

IF IOMP2 THEN 

H [ I , I - 3 ] : = 0.0; 



EMC ; (* LOOK *) 



PROCEDURE DOUBLEQR ; 
VAR T : INTEGER ; 

PROCEDURE COL. IC'D ; 
VAR I: INTEGER; 



+ABS ( R) ) ; 

Z) +ABS (il [ M + I / M t 1 ] ) ; 
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BEGIN (* COLMOD *) 

FOR I : =1 TO J DO 
BEGIN 

P := X* H [ I , K ] + Y* H [ I , K+ 1 ] ; 

IF NOTLAS THEN 
BEGIN 

P : = ? + Z*Z * : i [ I , K + 2 ] ; 

: I [ I , K + 2 ] := H [ I , K + 2 ] - P * R ; 

END; 

H [ I z X+ 1 ] := ii[I,K^I]-F*C; 

H[I,K] := H [ I , K ] - P ; 

END; 

ENT; (* COLMOD *) 

PROCEDURE RC'.IMCD ; 

VAR J: INTEGER; 

BEGIN (* ROvIMCD *) 

FOR J : = K TC G F 0 
E EG I lJ 

? := H [K,J] +G*H [K + l, J] ; 

IF NOTLAS THEN 
BEGIN 

P : = P + R * I : [ K + 2 , J ] ; 

H [ K + 2 , J ] := H [ X + 2 , J ] - P * Z Z ; 

END; 

H [ K + 1 , J ] := H [K+l, J]-F+Y; 

H[R,J] := • * [ K , J ] - F * X ; 

t . 

END; (* RONMOD *) 



PROCEDURE ACCTRANS; 

VAR I: INTEGER; 

BEGIN (* ACCTRANS *) 

FOR I : =LOW TO IGH DO 
BEGIN 

P : = X*Z [ I , X] +Y*S [ I , F + l] ; 

IF NCTLIG THEN 
BEGIN 

P : = P+ZZ*Z [ I , K+2] ; 

Z [ I , K+2] := Z [ I , K+2 ] -P*R; 

END ; 

Z [ I , K + 1 ] := Z [ I , K+l] -P*Q; 

Z [ I , X ] := Z [ I , K ] - P ; 

END; 

END; (* ACCTRANS *) 

EG I N ( * DOUBLEQR *) 

FOR K:=M TO NA DO 

B t_i O x 

* L fl — • 1_ • 1- / , I i ' j 

NOTLAS := KONA; 

IF KOI1 THEN 
BEGIN 

? := H [ K , K - I ] ; 



35 



tt tJ L 



Q := H [K+1,K-1] ; 

R : = 0.0; 

IF NCTLAS THEN 

R := H [ K+2 , K- 1 ] ; 

X := A3S(P) +ABS(Q) +ABS(R) ; 
IF X=0.C THEN 
FLAG1 :=' TRUE 

— i j ^C iU 

3 EG IH 

? := P/X; 

Q := Q/X; 

R := R/X; 

END ; 

END; 

IF FLAG1=FALSE THEN 
BEGIN 

S := P*P+Q*Q+R*F.; 

3 := SQRT(S); 

3 := SIGN (S , P) ; 

IF K= ! THEN 
uFGI M 

IF L<> 1 TIIET I 

H [ K , K- 1 ] := -H[K,K-1] 

SFD 

ELSE 

-J r rr T.r 11 . _ r * V • 

.1 l LK r L\~ JL \ . — ~ ^ / 

F : = P*fS ; 

X : = ?/ o ; 

:= C/S; 

Li ij . — l / t~J i 

Q := Q/P; 

R := R/P; 

P.CUMOD ; 

WRITE 

J : = MI NO ( EM , K + 3 ) ; 

GO L. 10D ; 

nrite 

AC CT RANG ; 

/RITE ( ’ . 1 ) ; 

END; 

END ; 

END; (* D0U3LEQR *) 

GIN (* HQRCPTl *) 

TSGTEXIT := FALSE; 

IF ITS>=100 THEN 
BEGIN 



I ERR 



:n; 



WRITELN ( ' THE ITERATION 
GRITEL.N ( ' REACHED IN 
SPACEBAR; 



= TRUi 



Li 



limit 

-HCR- 



OF ICO HAS BEEN 

) ; 



g i n 
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IF ( ITSOO ) AND (ITS MOD 10 = 0) 
BEGIN 

T := T+X; 

FOR I : = LCV7 TO Et: DC 
H 1 1 , 1 ] := H[I,I]-X; 

= ABF ( H [ EX , NA ] ) + ABF (H [ N 
= ''. 75 * 3 ;* 

= x' ; 

= - C . 3 7 5 * 3 * 3 ; 

E..D; 

ITS := ITS+1 ; 

LOCK; 

DCU3LEQR; 

END; 

END; (* HQRCPT1 *) 

PROCEDURE HCRCPT2; 

PROCEDURE F 0;DCC 2; 

VAR J: INTEGER; 

:ECIN(* RCNMCD2 *) 

FOR J : = I , A TC N DC 
BEGIN 

2Z := N[NA,J]; 

N[NA,J] := T.*ZZ + P*H [EN, J] ; 
H[Ei:,J] := r*H [EN, J] -P*EZ ; 



END; (* ROVJIR 



PROCEDURE CCLNCD2; 
VAR I: INTEGER; 

BEGIN (* COLMCD2 *) 
FOR I 



; ( 



: =1 


TC 


D7 7 T 
i-J L 4 


DO 








-IN 














n . 


= H [ 


I,N 


-] ; 








[ I , 


NA] 


; = 


r\ z 


n 

u 


+ P 


* N [ I , EN ] ; 


[I, 


EN] 


: = 


0* H 


r 

L 


I , 


EN] -P*ZZ ; 


COL 


^ ,r ~ -s o 

> iC U « 


*) 











PRCCEDURE ACCTRAN32 ; 

VAR I: INTEGER; 

i EG IN (* ACC TRANS 2 *) 

FOR I : = LC N TO IGH DC 
BEGIN 

ZZ := Z [ I , I! A ] ; 

7 [ I , N A ] := Q* Z Z + P* Z [ I , EN ] ; 

Z [ I , EN ] := C* Z [ I , EN ] - P* Z Z ; 



E 




^ * 



) 



(* HQROPT2 *) 
(Y-Xl/2.0; 



THEN 
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Q := P* P+W ; 



ZZ := SORT (ADS (Q) ) 


H[EN,EN] := 


X + T ; 


X := H f E X , E N 


i • 

i / 


H [NA, XA] := 


Y + T; 


IF Q < C . 0 THE 


LI 


3 EG IX 




••E [NA ] : 


= x+P; 


'-[NX] : 


= X + P ; 


.71 [XA] : 


— LJ f 


■ 71 [LX] : 


= - ZZ ; 


END 




ELSE 




BEGIN 





ZZ := P + 3 IGM ( Z Z , P ) ; 
NR[NA] := X + ZZ; 







[El ] 


; = 


NR [ NA 


] ; 


I 


P 


ZZOO . 


0 THEN 








•7R [ E 


N 1 


• - V _ ‘ 


/ZZ; 


< 


I 


[NA] 


: = 


0.0; 




. i 


T 


[EX] 


: = 


0.0; 








Cl! 

II 


t NA ] ; 




o 

i\ 




: = 3 


m 

A 


(X*X+Z 


Z*ZZ) ; 


p 




• - V 


/R; 










. _ r? 


Z/R 


r 




r.v 

L \ 


r\ 


:?mgd 


O . 

“ f 






r 


R 


T m t- / 

± i IL ^ 


1 1 


) ; 






r\ 


L.ICD 


Z } 






. » 


iv 


ITS ( 


t i 


) ; 




A 




CTRA 


NS 3 


r 




* 


p 


ITZ( 


1 1 


) ; 





-xo ; 



i Zjl t . — Jj IN 1 1 4 f 

IF EN> = LO.J THEN 
BEG IX 

ITS := 0; 



i—i i i !_/ 

ELSE 

FLAGS := TRUE; 

EXD; (* HQRGPT2 *) 

PROCEDURE EIGENVECTOR; 
VAR EX , I , J : INTEGER; 

PROCEDURE LCCP700; 
VAR I , I I , J : INTEGER; 

TTP : REAL; 

BEGIN (* LCCP70C *) 
WRITE 

FOR 1 1 : =1 TO N A DC 
BEGIN 

I := E X -II; 

,1 := H [ I , I ] - P ; 



38 



R : = H [ I , Eli ] ; 

IF M< = NA THEN 

FOR J:=M TO NA DC 
BEGIN 

IF ( H [ I , J ] =0 . 0 ) OR ( H [ J , EN ] =0 . C ) 
TTF := C.C 



THEN 



TTF := H[I,J]*I 
R : = A - T T P ; 

Li 4 f 

if '.:i [ i ] <c . o the:, 1 

DEG I LI 

7 7 • = ;7 • 

— ^ " ' * r 

r* . — 






: = r. 



END 
ELSE 
3 EG I M 
;i 



T . 



IF ,71 [ I ] =C . 0 THEN 
BEGIN 



. _ ' 7 . 



l T ■ ’ = Q Q 'THEN 



- ^ ghep*::or:i; 



II [ I , EN ] := -R/T; 



ELSE 

LOGIN 



i: [ I , i + l] ; 

H [ I +1 , i ] ; 

( NR [ I ] -P) * (NR [ I ] -P) +.71 f I ] * 
(X*£-ZZ*R)/C; 

H[I,EN] := r; 

IF AO S ( X ) > A3 3 ( Z Z ) THEN 
H [ I + 1 , EN ] := ( - R-v7*T ) /X 

ELBE 

H [I+l, EM] := (_s-Y*T)/ZZ; 



71 [I] ; 



’::rs • 



END; (* LC0P70C *) 



PROCEDURE LCCP79C ; 
VAR 1 1 , J : INTEGER; 



PROCEDURE LOCPWORK; 



GIN 


(* LOO PNC 


EK 


*) 


X : 


= H [I, I+l] 


t 




Y : 


= II [ I + 1 , I ] 


r 




VT\ 


:= (NE [ I ] - 


F) * 


(NR [ I ] -P) +V7I [ I ] *;7I [ I ] 


f - T 

V JL 


II 

J 


? ) * 


2 . c * c ; 


IF 


(VR=0.G) A 




(VI=C.O) THEN 



\ V U I ^ j i - «. . — / \ V L — • yj J I ULl, 

VR := MACHZP*NORN* (A3S (W) +AC3 (C) + 
ADS (X) +AEC (Y) + ABE (ZZ) ) ; 



G C M D I V ( X * R. - Z Z * R A + Q * 3 A , X * S - Z Z * 3 A - Q * R A 
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VR, VI , ZR, ZI ) ; 

H[I,NA] := ZR; 

H[I,EN] := ZI; 

IF A 3 S (X) > (A 3 S (ZZ) +ABS (Q) ) THEN 
BEGIN 

i![I + I,MA] : 

H [ I + I , EN ] : 



T V T 

‘-.JVJXU 

CCilDIV (-R-Y*H [ I , N A ] ,-S-Y*H[I, EN] , Z Z , Q , Z P , Z I ) ; 
H[I+l,NAj : = ZR; 

H [ I + 1 , EN ] := ZI; 

END; 

ID; (* LOOF'JORK *) 



= ( - RA- VI* H [ I , NA] +C*H [ I , EN ] ) /X ; 

=* (-SA-VI*H [ I , EN] -t*il [ I , NA] ) /X ; 



BEGIN (* LCOP790 *) 

FCR 1 1 : = 1 TO Ell; 1 2 DO 
BEGIN 

••RITE ( ' . ' ) ; 

t • _ v; i t t . 

-L • — Ll.l“ 1 1 ; 

.1 := H [ I , I ] - P ; 

RA : = 3.0; 

SA := II [ I , EN] ; 

FCR J : = M TC NA DO 
BEGIN 

RA : = F A + H [ I , J ] * F [ J , NA ] ; 

3A := SA+H [ I , J ] * H [ J , Eli ] ; 

END; 

IF . J I [ I ] <0.0 THEN 
BEGIN 

r-» f» . r T . 

Li • — r V f 

R : = RA ; 

G := SA; 

END 
ELBE 
BEGIN 
N : = I ; 

IF J I [ I ] = 0 . 0 THEN 
BEGIN 

CONE I V ( - R A , - 3 A , N , C , Z K , Z I ) ; 
H [ I , NA ] := ZR; 

H[I,E1I] := ZI; 

END 

ELSE 

LCOPWCRK ; 

END; 

END; 

END ; ( * LOOP790 *) 



PROCEDURE LCOPG30; 

V \R I , J , K : INTEGER; 
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BEGIN ( * LOOPS S 0 *) 

FOR J:=N DOWN TO LOW DO 
BEGIN 
WRITE 

M := MI NO ( J , IGH ) ; 

FOR I : =LOW TO I GW DO 



7 1 ' ]? T • — ► ! : 

L ^ L\ f\ * — ^10 J l .1 lJ w 

ZZ : = ZZ + Z [ I , K] *H [K, J] ; 

Z [I, J] := ZZ; 

END; 

END; 

END ; ( * LOOPS 30 *) 



- - ( ri [ E N , EN ] - P ) /H [ E , IA] ; 



PROCEDURE PRELOOP ; 

BEGIN (* PRELCOP *) 

' * — • 

IF A33 (H [EM , NA] ) > A3S ( !! [HA , EN ] ) T'lEN 
BEGIN 

N[N0 ,!?:.] 
i! [ NA , EN ] 

lli D 

ELSE 

BEGIN 

BOND I V ( 0 . 0 , - H [ NA , EN] ,11 [ NA , NA ] - P , Q , Z R , Z I ) ; 
H L l\ A f ] l = N R ; 

H [:ia,e::] : = ZI ; 

END; 

u r TT • * T l 1 . - n • 

II [ EN , EN ] := 1.0; 

t? m > * ^ ^ *1. 

Ui 4 t . W • An*' 1 ; 

END; (* PRELOCP *) 



GUT ( * EIGENVECTOR *) 

NORM : = 0.0; 

K : = 1 ; 

FOR I : = 1 TO N DO 
BEGIN 

FOR J : =K TO N DO 

NORM : = H0EM+A3S ( H [ I , J ] ) ; 
K : = I ; 

END ; 

IF NORMOO.C THEN 
3 EG IN 

FOR EN : =M DOWNTC 1 DO 
B EG I N 

P := W 7 R [ EN ] ; 

0 := WI[EN]; 

NA := EN-1; 

IF Q<=C.O THEN 

.o vj 1 a 

IF Q = 0 . 0 THEN 
EEC IN 

'4 . — p M . 

1 4 •“ L i, « r 
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H [ EN , EN ] := l.C; 

IF UAOO THEN 
LOOP700 ; 

END 

ELSE 



.'I 

PRELGCP ; 

-l. s ill — .; 

L..OP730 ; 

END; 

END J 

tt» ' T r> - 

FOE I : = 1 TO M DC 

IF (KLCW) CF. ( I > IGH) THEN 
FCR J : = I TO N DO 

Z [ I , J ] : = H [ I , J ] J 

LCCPS80; 

END; 

■?Mn • /* "7 t p \tt / ^ m p r-' x \ 

<-• l * 4 -/ / { ui ) 



IIN{* HC.R2 *) 



!R I : = i TC N DC 
BEGIN 

IF (KLCL) OR ( I > IGH ) THEN 
BEGIN 

•R [ I ] : = H [ I , I ] ; 

..’I [I] t = 0.0; 

END; 



:= IGH; 
: = 0 . C ; 



"m s = r ^ 



THEN 



! T ? n T ^ t 
t VJ 1 IT* 



= 0; 



MA := EN-I; 

ET 12 := IiA- 1 ; 

Tr . _ ir'7'r^p. 

RPSAT 

FLAG 4 := FALSE; 

L := ZN + .1 ; 

REPEAT 

L := L- I ; 

IF L = LO’N THEN 
FLAG4 := TRUE 



IF ABS (H [ L , L- 1 ] ) <=.'1ACHEP* (ADS (H [ L- 1 , 
) THEN 

FLAG 4 := TRUE; 

UNTIL FLAG4-TRUE; 



r . _ »J Tr'T T? ' 7 1 

‘ • - il l L -i , ZjL'i J 

:f loen then 

j. E G I N 
Y : 



H [ NA , NT. ] ; 

H [ EN , NA] * !! [ NA , EN 



j '/ 



L-l] ) +ABS (H [L,L] 
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u rrc 



IF LONA THEN 
BEGIN 

HCR0PT1 ; 

IF TESTEXIT THEN 
EXIT; 

END 

ELSE 

HCECPT2 ; 



1 L S E 
BEGIN 

H [ EM , EN ] := X + T; 

NR [ EM ] : = H[EN,SN] ; 

wl [EN] := O.C; 

EM := NA; 

IF EH>=LOW THEN 



i :i 

ITS 
NA : = 
ENN.2 
END 



1 - 

JC. J- f 



' _ 1 • 
— t 



? r A n *3 • — ff D r 7 n . 

— vj • X j. v ^ i_j j 



END; 

' I L FLAGS ; 



i i L 
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(* VERSION 0206 *) 

MODULE OVERLAYS; 

(* MODULE EIGEN3AL *) 

TYPE DCMAINl = 1..2C; 

MATRIX = AP.r '.Y [D-: LAI Ml , DOMAIN:. ] CP DEAL; 

LIST I = ARRAY [DOMAIN 1] OP INTEGER; 

EISTR = ARRAY [DOMAIN 1] CP REAL; 

PROCEDURE EL LUES (MM , N , LO'M, I GH : INTEGER ; VAR A:MATRIX;VAF IiIT:LISTI) ; 
VAR I , J , KPl , LA , M , MMl , MPl : INTEGER; 

X , Y : REAL; 

BEGIN (* 2LI1HES *) 

LA := I GH- 1 ; 

DPI := LC.7+1; 

IF L A > = K P I THEN’ 

ELGIN 



MMl := r-I-1; 

X : = 0.0; 

I : = 

FOR J : = 1 TO I G : I DO 
ELGIN 

IF AD 3 ( A [ J , MM 1 ] ) > AB G ( X ) THEN 
BEGIN 7 

■ •’LITE ( 1 . ' ) ; 

X := A [ J , MMl ] ; 
i. . U / 

2ND ; 

END ; 

I NT [M] := I; 

IF 1.0 l THEN 
BEGIN 

FOR TO M DO 

R EG I N 

Y : = A [ I , J ] ; 

A [ I / J ] := A [ M , J ] ; 

A [ M , J ] := Y; 

END; 

FOR J : =1 TO IGH DO 
BEGIN 

Y := A [ J , I ] ; 

A[J, I] : = A [ J , M ] ; 

A [ J / M ] : = Y ; 

END; 



IF XOC.C THEN 

O T-.Q I I 

MPl := 1+1; 

PCR I : =MP1 TO IGN DO 
D E G I N 

Y := A [I, MML]; 
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LU 



r-\ 

U 



L.N 1) ; 

END; 

SLMHES * 



IF YOO.O THEN 
BEGIN 



TRITE ( ' . ' 


) ; 




Y := Y/X; 

A [ I , MM 1 ] 


: = 


Y; 


FOR J : =M 


T 3 


N DC 


■A [ I , J ] 


: = 


A[ I , J] - Y * .* [M, J] ; 


FC R J : = l 


TO 


IGH DC 


A i J , .’•! j 


: = 


\[J ,-A] +Y*.. [ J , I ] ; 



END; 



OCEDURE ELTRAN (N 1 , M , LON, IGH : INTEGER ; A: MATRIX; 

JN : L 1ST I ; VAR Z : MATRIX ) ; 

R I, IF,MF1: INTEGER; 



G I N( ' 



RAN *) 



DO 

a , 



FOR I : = I TC i! DO 
BEGIN 

FOR J : =1 TO I 

z [ I / J ] : = 

Z[I,I] : = 1 . 0 ; 

Ij } 

KL := IGH-LCN-1; 

IF KL> = 1 THEN 
S E G I N 

FCi< HM:=1 TC XL DO 

£3 Jj ‘o I i. « 

M? := IGH-MM ; 
MPl := . M P + 1 ; 



FOR I :=:!?! TC 



T -I! DO 



Z[I,/.F] := A [ I , MP- 1 ] ; 

I : = JN [Mi?] ; 

•TRITE 

IF 10:1? THEN 
BEGIN 

FOR J : =MF TC I Gil DO 
BEGIN 

Z [ M P , J ] : = Z. [ I , J ] ; 

Z [ I , J ] : = 0 . C ; 

END; 

Z [ I , IIP ] := i.C; 

END; 



Nil 



D; (* ELTRAN *) 



OCEDUFE BAEBAX (NM, N, LON, IGH : INTEGER ; SCALE : LISTR ; 

Mi : INTEGER; VAR Z MATRIX) ; 
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INTEGER 



VAR I , J , K , 1 1 : 

S: REAL; 

3 EG I M ( * 3AL3AK *) 

if igholcg the;: 

3 EG II! 

FOR I : = LC . rc IG” 'l>C 
3ZC-IG 

J := SCALE [I]; 

F SR J : = 1 TO '1 DC 

Z[I,J] : = Z [ I , J ] * S ; 

END; 

EGD ; 

FUR 1 1 : =1 TC G DO 

n p t »t 

oL'Jiii 

I : = II; 

IF (KLC-vJ) CPs ( I > IGH ) THEN 
3 CC I U 

IF KLCN THEN 
I := LC»;-II? 

:= tru-;c(s:alc[i] } ; 

IF KOI THEN 

■: eg i 

;FITE 

FOR J : = 1 TO ;• DO 
IF GIG 



2 [ I , J ] 
2 [ K , J ] 



[ i , J ] ; 

• = -j [ K ^ j ] ; 






END; 

EHD;(* 3AL8AK *) 



PROCEDURE BALANCE ( U 1 , M : IUTECER; VAR o : MATRIX ; VAR 

VAR SCALE : LI 3TR) ; 

CONST RADIX = 2.0; 



r - - T ^ r J • T' Tr P' 
IjO rv f 1 A . i o 1 . 



VAR FLAGA , FLAGD , FLAGG , FLAGD , FLAGI , FLAGS , GCCCGV : 
I , I EXC , J , K , L , • I : I GTE GER ; 

32,C,F,G,R,3: REAL; 



B00LEA: 



PROCEDURE RC Go E ARCH ; 
DEGIG (* RC'GGEARCr! *) 
J := L+l; 

FLAGA := TRUE; 
REPEAT 

J • \J ~ _L / 



^ L . . j 4 . *— i' j'. kjo Tj ; 

REPEAT 

I : = I + 1 ; 

IF IOJ THEN 
8 E G I G 
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IF A [ J , I ] <>0 . 0 THEN 
FLAG2 := TRUE; 

END ; 

UNTIL ( I =L ) OR ( FLAG2=TRUE ) ; 
IF NOT FLAG2 THEN 
BEGIN 



I EX( 



1 1 
' , Mr> , 



J IA TIL 

END; (* 



(J=l) OR (FLAGA=7AL3E) 
RCNSEARCII *) 



PROCEDURE C 0 L U ' -1 N 3 E A R C H ; 
BEGIN (* COLUMNS E ARCH *) 
J := K- 1 ; 

-LnCB * “ TRU~ ; 

REPEAT 

J := J+l; 

I : = K - 1 ; 



FL 

RE 




FALf 



i . - 1 TJ. , 

IF IOJ THEN 
BEGIN 

IF \ [ I , J ] < > G . C THEN 
FLAC-1 : = TRUE; 

END; 

UNTIL ( I =L ) OR (FLAG1=TEUE) ; 
IF N'CT FLAG1 THEN 
BEGIN 

N := N; 

IEXC := 2; 

FLAGS := FALSE; 

END ; 

UNTIL ( J = L) CR. (FLAGB = FALSE) ; 
END; (* COL U\ INSEARCH *) 



PROCEDURE NCR. ‘.REDUCTION ; 

VAR I , J : INTEGER; 

BEGIN (* UCRHR EDUCTION *) 

FOR I : =K TC L DC 
BEGIN 

C := 0.0; 

R : = 0.0; 

FOR J : =K TO L DC 
BEGIN 

IF JO I THEN 
BEGIN 

C := C + ABS ( A [ J , I ] ) ; 
R := R + AB.U (A [I , J] ) ; 
END; 
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G := F/RADIX; 
WRITE 
F : = 1.0; 

S ! = C + R ; 
WHILE C<G DC 

3 r.^ [ m 




G := R* RAD IX; 

WHILE C> = G DO 
3EGI tl 

F : = F/RADIX; 

C := C/32; 

END; 

IF ( (C + R)/F) <(0.95*S) THEN 
BEGIN 

-j * - 1 • 0/ c ; 

SCALE [ I ] : = SALE [ I ] * F ; 

•’/> . mT'f’r , 

. 4. t J • i l \ J L ^ 

TOR J:=X TO N EC 

A [ I , J ] : = A [ I , J ] * G ; 

FOR J : = 1 TO L CD 

A[J,I] := A [ J , I ] * F ; 



END; (* NOR.' DEDUCTION *) 



BEGIN (* BALANCE *) 
BE := RADIX* RAD IX; 

X : = 1 ; 

r • = • 

* ' * ! 

y? ^ • — i? a r o’ r . 

L u ri O O * — l nL». _j f 

WHILE NOT FLAGD DO 



BEGIN 
P L oD 



:= TRUE; 
:= FALSE; 



-T • 

• * / 



IP PLACA fHEII 
COLUMNS CARCH; 

IF HOT FLAGD THEN 
3 EG IN 
REPEAT 

SCALE [M] := J; 

IF JO.l THEN 
BEGIN 

POP I : =1 TO L DO 
BEGIN 

? : = A [ I , J ] ; 

A [ I , J ] := A [ I , : 1 ] ; 

Mi/'i] : = - ; 

END ; 

FOR I : =K TO N DO 



IP/^T ' ’ 
) L G 1 il 



= A [ J , I ] ; 
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A [ J , I ] : = 

A[M,I] := 
END; 

END ; 

IF I EXC= 1 THEN 
FLAGC := TRUE 
ELSE 

t- ->\ . _ T7 A f r’ rp . 

L Li r - J -w • “ - rn i_> *.> ^_i f 

I? NOT FLAGC THEN 
DEC IN 

K := K+i; 
COLUMNSEARCH; 
END; 

.71' I TE 

UNTIL FLAGC OR FLAG 
IF NOT FLAGB THEN 
BEGIN 

IF L<> 1 THEM 
BEGIN 

FLAGC : = FA 
L : = L-l; 

/ 

END; 

END; 

END ; 

IF FLAGB THEN 

^ T ' ■ 

77 o r> t . — r - .no r o o 

i- VD v 1 • i \ x v> Ls 

SCALE [I] := 1 . C ; 

REPEAT 

• ; G PM REDUCTION; 

UNTIL NOT NOCCHV; 
iii N D j 

LOW := X; 

I OH := L; 

ELD; (* 3AL.VJCE *) 



XCDEND. 






A[M,I] ; 
F; 



3; 
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(* VERSION 0235 *) 
(* MODULE DETERM *) 

MODULE OVERLAY 10 ; 



. _ -v 



CM'. I Ml = !. 



’ 1 ^ M A 

^ r ^ 



T ? 1 p ' 
. i. j 



pf^mj _ a ^ i • " 77 1 

Lj ^ o i .vr. . [ . •. j. ■» ^ ^ i. 

LI3TR = ARRAY [DCMAINl] OF REAL; 



VAR IOMCIU: EXTERNAL I iTEGER; 

AO I FLAG : EXTERNAL 2 CO LEAN ; 

IOTA: EXT JAL MATRIX; 

'■* c r> tt •• "> • t? "/ m ^ p ' ■ a r n r 1 ' r . 
j C i • L, A i ill I v - A Li r\ j_. c \ L f 

EXTERNAL PROCEDURE CLEARSCREEN; 
EXTERNAL PROCEDURE SPACEBAR ; 



EXTERNAL FUNCT 
EXTERNAL [2] F 



iCED jRI 



DOC LEAN; 

MAFOCCPY (IINUM: INTEGER) 



EXTERNAL [1] PROCEDURE RLUD ( ND , N : I NTEGER ; VAR KER: 

MATRIX; 

VAR JN 



rv • r.ix 



p p C'. r* ^ y p p p m M'vET * 

rj > p'. -» i-» . ^ppp . 

r- p p-i •" f r> • r t -p • 

l - U Li • LJLI-DX^. f 

"I T ' T • T T C -7 1 T • 

L JL * « • i_l _ i J. f 

3DET : REAL; 

I , J , PDEC , NT : INTEGER; 

CMC ICE, SELECTION : CHAR; 

r->/-~',r>T-p # r- r ('-> m p . 
j. ^ l. • L i v i r. / 

r ■ : r liztr; 

i : iijtege::; 

p ! T ■> TN • 

- “ ■ ^ w £ * A / 

OOrnf . ■ ’n’n T ‘P • 

*— i.\ <— > > j — j • vJ X A 1 Nv.) ^ 

D * . T ' Tf-l'T'/"* nn . 

r^. • 1 h i liVw n, K f 



l . U .. 1 -Mi L\ 



I? T N 



PROCEDURE F.DET ( SND , EN : INTEGER; VAR OREL: INTEGER ; VAR 

VAR 3A: MATRIX) ; 






VAR J: INTEGER; 

3 J I : LI ST I; 
33CALE : LIETR; 



.CUii.i J) f 



CM: CMC R ; 



BEN, IN {* RENT *) 

SEER := 0; 

RLUD (SNE , Sti , SKER , E 



A , E J N , NSC ALI 
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DET := S JM [SN] ; 

FOR J : =1 TC SC DO 
DET := DET*SA[J,J] 
CMD ; (* ROOT *) 

2ZGIC (* TTFDE7 *) 

r* r - ^ . 

•• — -• - ' - • 4- -> «/ 

:~i . _ jQ"rs r j . 

r > — 1 T . _ 1 .-*o > rr •^i"' 

L' » v . — . . , . i ^ ^ 

FOR J : = 1 TO GT DC 



1 r ' r t 1 
t V Lj ^ -i. f ^ \ 



OTA [ I , J] ; 



i \ . 

t t 



CLSARCCREZII; 

HEITELH (' PLEASE JA I T 
REST ( MT , II T , TKER f SDET , TALC ) ; 

r TT 7\ TO C r ^ TR * ~ . 
v-iJ Li Ui, f 

IF TKER = C the:: 



'7P pryi . 

— jD L i ^ 



14 Ji J 

^ r>r\ T'm 
OoO i-« 1 

./3ITEL:; ; 

: .;r.:TEL: T ( 1 t:ie dete :? = 1 ,sczt) ; 

n p \c PBTi r 

7 -p T rrT -1 - « - . 

4 - . 1 X w LJ M r 

'TDTn;' / I v ''0 r 7 T ■ a ' ’ n a :jr oo ^ ^ o/ *2 ■» r / m *\ * 

* t\ x j. l- ^ X. ' j U . J n : . x 1 1 1 1 r. . \ iv u. x. . ^ / l . ; r 

IF YES THE'! 

HARDCOPY ( 5) ; 



e::d; 

FOOD ; ( * 



100 END . 
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(* VERS 10 'I 30 0 *) 
(* AX=B MODULE *) 
MODULE OVERLAY 20; 



TYPE DOM A I Ml 
MATE I ' 
LI JTI 

f T r-! 



= 1 . . 20 ; 

ARRAY [DCMAIUl ,DCHAIi:i] 
T \Y [ DC MA I’M L ] OR IMTEG 
OR PERL; 



p ^ ' r • 

: i_ . r 



v r ; 



r . T v 1 1 



ict; 

AC II 



EXT 



Lj i\ x i-t 



' - r r n *— -> • 



, ICTB , XX : EXT E REAL MATRIX; 

'LAG, IODRLAG: EXTERNAL BOOLEAM; 
1:1, I CM: EXT E DUAL INTEGER; 

REAL FUMCTICM YES: BCOLEAU; 

REAL FRCCEDURE SPACEBAR; 

CLEARSCREEr ; 

UR E RLUD 



( 'T 



"* “R • T/A TO • T : T ^ ^ P • 1 

. f v »l\ . i. — _ L f 



;xt: 



'ATE IX; 

VAF J:;:LISTI ; VAR SCALE: LISTF. 
FIIAL [1] FRCCZDUFE RFRS (VC, R: INTEGER; VAR KZR : I MTEGER ; ' 

VAR JU: LISTI ; VAR X : LI STR) ; 



’ERUAL [-] PROCEDURE HARDCOPY ( HMUH : I MTEGER) ; 



r-n r* -a < ' • 

i 1 ^ »R A_ f 



AR QUITFLA3: BCCLEAM; 

I , J,PSE3: IMTEGEF; 

CHOI IF, SELECT I CM: CHAR; 



T P* -p .If* * T - T >. t n " r> -3 . 

L Lj L i , k 1 j. j l* ± m X . . ± u'JUl f 

CM: CHAT ; 

SREEL: STR I MG; 

?4: I MTEGER ; 

rp T, RO • T>’rnRC^n . 
x_.j~.~v. 1 A x L O .j A t - 

i u.. 'w j. .* Li i j • Lj .l. i i \ ^ 

TIM: LISTI; 

TIM IT: I MTEGER; 

TALU,TB: MATRIX; 



PRCCEDUr E SAXB (SMD, SM,, 



M : I MTEGER ; VAR G A : MATRIX ; VAR 
VAR GIMIT : IMTECER; VAR S JM : LISTI ; VAR 
AAR I , J , M M : I MTEGER; 

3B , SC. LE : LICIT; 

3 EG I C (* SAXE * ) 






1 i 



:mit=o them 
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11 Cl 



RLCJD ( SND , SN , SKER , SA , S JN , SCALE) ; 

IF S JIT [ StI ] = 0 THEM 
BEGIN 

SKER := 3; 

WEITELN; 

XETTELN ( ' IN SAXE , L’J DEG OVIPOSIT I 
NEITELN (' YIELDED- A SINGULA? U . 

if.iteln ( ' solution deed not exist 

■.'o' n rc\ r ' • 

x f 

EXIT; 

END; 

"■I? r*.i / - A rr.r - 
1 t 1 \ - U i 11 L. H 

c IER := C; 

EXIT; 

END; 

EIR J : = 1 TO 3 VI DO 

nri*' t », * 

E i'j .i x * 

fcf i : = l tc :::j cg 

DD f T 1 • - f T T1. 

L -L J • >- iJ t ^ f J J / 

n TT 1 ”" ’ / '" * T r" • ■» ■- r* r-1 — > ’A c ' r ' ’ ~ d \ • 

- o \ -D n u / w i ^ w i.\ i j :\ r * - -i i Oa f .-j ) { 

FOR I : = 1 TC JN CO 
3 X [ I , J ] := ED [I] ; 



END; 0 



) 



GIN ( * TTEAXD *) 
CLZAR3CR EE A ; 



- t 






TIN IT := r ; 

FOR I : =1 TO .IT DO 
FCR J : = 1 TC NT DC 

TALU [ I / J ] := IOTA [I, J] ; 

FCR J : =1 TO IT DC 
FOR I : =1 TO RT DC 

CLEAPSCREEN; 

.JRITELN (’ PLEASE XAIT ’); 

3 AXE ( NT , NT , NT ; TALU , T3 , XX ,T IMIT , TIN , TKZR 
CLE.ARSCr DEN ; 

IF TKDP=0 THEN 

FOR J : =1 TC IT DC 
3 EG IN 

’..'El TELK ( ' THE SOLUTION FCR CCLUM 
FCR I : =1 TO NT DC 
BEGIN 

SUITE ('['); 

NDITE (XX [ I , J] ) ; 

NRITELN ( ' ] ' ) ; 

END; 

EP7.GD3AI ; 

NRITELN; 
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* u 



mi JD ; 

WRITELN ; 

WRITE (' DC YOU WANT A HARDCOPY? Y/N ' ) ; 
IF YES THEN 

JARDCOPY ( 2) ; 



ICDELID . 
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(* VERSION 1283 *) 

(* POLYROOT MODULE *) 

MODULE CYSKLAY2 1 ; 



CCNJT MAXDEGPl 
RAXDEC = 



= 21 ; 

?0 : 



:EE DO: :A I Ml = 1..20; 

DGlIRPI = 1 . . MAXDEGPl ; 

MATRIX = ARRAY [DCMAIM1 , DOMAIK1] OF REAL; 
LIETI = ARRAY [COMA I XI] OF INTEGER; 

LIETR = ARRAY [DOMAIU1] OF REAL; 

LISTRF1 = ARRAY [COMFPl] OF REAL; 

EXTERNAL PROCEDURE CLEARSCREEN; 

EXTERNAL PROCEDURE CLEAR IT ( I : INTEGER) ; 

EXTERNAL PROCEDURE IN TREAD (VAR K: INTEGER) ; 

EXTERNAL PROCEDURE SPACEBAR; 



V m p r; r. - - t [''?"] ; 

^ A 1 IS illDJ MI U [ S J £ :R uci. U i\ I 



r 7 . t ’ T rr' • TO . T ? " T. 0 ~ P 07 • 

jL 4 r' v_/ J v « i\ . J . i- 1 < i ~ . j ~ . , V -A I igUu- . 



EXTERNAL [3] PROCEDURE HCR 2 ( MM , H , LOOJ , I Gil : I NT 

EE, Ml : 



VAR M : 



, i o i 



T? • 3 * ^ p 7 . M ^ f? E T V • r TEO 

1 i \ : \ ^ • mm 1 i 1 f V / i ,l uiv 



EXTERNAL [6] PROCEDURE BALANCE ( NM , N : INTEGER; VAR A: MAURI 

IGH : 

INTEGER; VAR SCALE: LIST 1 ? 

EXTERNAL [5] PROCEDURE ELMHES ( 'R ! , U , LCN , IG:i : I NT EGER ; V'R 

VU R TNT: 

LI ST I ) ; 

EXTERNAL [6] PROCEDURE ELTRAM ( NM , N , LON , IGH : INTEGER ; VAR 

JN : LIST I; 

VAR: Z: MATRIX); 

PROCEDURE TTRPQR ; 

VAR QUITFLAG: BOOLEAN; 

I , J , PDEG : INTEGER; 

CHOICE, SELECTION : CHAR; 

PRCCEF: LISTRP1; 

TMI,T’;R: LICTR; 

T I ERR , NDEG : I NT EGER ; 

C " • p M A O • 

~ ‘ • l ' / 

oRi-ijL: GTR i ij>j j 

P4: INTEGER; 

TABGERR , TRELERR : LISTR; 
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TKLUST : LISTI; 

'll , TKER : INTEGER; 



PROCEDURE P 3 '.ID 2 (Ri: : INTEGER ; COIF : L I STEF 1 ; P. ■!?, , RE ! : LIGTR; 

VAR A S'S ERR , RELERR : L I STR ; VAR K LUST : L I ST I ; 
VAR XER: INTEGER) ; 

VAR I , J , J F 1 , J F , K , M , :.IM 1 , N P 1 : INTEGER; 

E , CERT D I ST , EM AC , MAG , OLDERE , P , PCNER , F. , RAT , 

SVR , SVI , Jl’CERT , VP. , VI , VT , X. I AG , XI , XR : REAL ; 

SHRUNK: BOOLEAN ; 

FUNCTION CPASS (XR, YR: REAL) : REAL; 

BEGIN (* CPA3S *) 

CPA33 := SQRT ( XR* XR+YR* YR) ; 

2ND; (* 



AMAXI (A ,3: REAL) 
3 EC IN (* AMAXI *) 

IF A<3 THEN 
AMAXI := 3 
ELSE 






A ;A X 1 



FUNCTION SINGLE (VALUE: REAL) : REAL; 

ELGIN (* SINGLE *) 

SINGLE := VALUE; 

ELD ; (* SINGLE *) 

FUNCTION DOUBLE (VALUE : REAL) : REAL; 
BEGIN (* DOUBLE *) 

DOUBLE := VALUE; 

END; (* DOUBLE *) 



PROCEDURE SECNALF; 

VAR J , K : INTEGER; 

BEGIN (* SECNALF *) 

FOR J : = I TO RN DO 
DEO IN 
TRITE ( 

X L U 3 T [ J ] := 1; 

XMAG : = GRABS ( RWR [J] ,ETI [J] ) ; 
EMAG := AB3ERR [ J ] ; 

IF EMAG=C . 0 THEN 

R : = 0.0 

ELSE 

ELGIN 

IF X: LAC = 0 . 0 THEN 

E : = -1.0 



EL Oi 
K 



/"M’' ^ t 

— vJ/ /V. 1/iLj f 
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END ; 

RELERR [ J ] := R; 

END ; 

NH1 := RN-I; 

FOR J : =1 TO NM1 DO 
BEGIN 

;kite 

J?1 := J+l; 

FOR X:=JPi TO ?.N DC 

oEjit* 

dist := .:?ab5(Rw’r[j]-r;;r[:<] ,i 

IF DIGT<= (ABSERR [ J] +A3SERR [K] 
BEGIN 

KLUST [ J ] := XLUST [ J ] + 1 ; 

KLUST [K] := KLUST [K]+l; 

END; 

END; 

END ; 



i.r. a 

:d; 



1ECHAL 



) 



:GIN (* REND 2 *) 

NRITELN; 

IF RN< 1 THEN 
BEGIN 

r -j ’a * = — ; 

NRITDLN (' N (DEGREE) MUST BE 
SPACEBAR; 

EXIT ; 

rr*’n . 

L ji< J / 

NP1 := RN+1 ; 

POWER := 1.0/RN; 

P := '33 (COEF ( 1 ] ) ; 

IF P=0 . 0 THEN 
BEGIN 

KER := 2; 

NRI TELN' ( ' LEADING COEFF I G I ENT 
SPACEBAR; 

EXIT; 

END ; 



v _ "I I 



) ; 



RA 1 ! 



: = 


AB 5 


(COEF [ 


2 JP 1 ] )/P? 


: = 


RAT 


* EXP ( ( 


- 45 . 0 ) * L 


JR: 


= 1 


TO RN 


DO 


i G I N 








WRITE ( 


' ' ) 


/ 


XR 


: = 


DOUBLE 


(R/JR [ JR] 


XI 


: = 


DOUBLE 


(RNI [JR] 


VR 


: = 


DOUBLE 


( 0 . 0 ) ; 


VI 


: = 


DOUBLE 


( 0 . 0 ); 


TOR 


J : 


= 1 TO 


NPl DO 



:gim 

WRITE ( 1 
VT : 

VI : 



ili L 



X R* VK- X I * V I 4- DOUBLE (CC; 
X R * V I + X I * VR ; 



[: 



Nil [J] -RNI [K] ) ; 
) THEN 



vC . ’ ) ; 
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a Oi 



VR : = VT ; 

END; 

3VR := SINGLE (VR); 

SVI := SINGLE (VI); 

NAG := CPA33 (SVR,SVI) ; 

3 : = A ■ i AX 1 ( RAT , ‘ 1AG/F ) ; 

RELERR [JR] := 1 ;* 

AES ERR [JR] := EX? '?0..'ER*LN (3) ) ; 

END; 

JNPUNX := FALSE; 

PE FEAT 

/rite 

SHRUNK := FALSE; 

FCR J : =1 TO PM DO 
3 EG IN 

IF A3SERR [ J] <>C . 0 THEN 
BEGIN 

:i : = RN; 

FOR ” : =1 TO EM DC 
DEG IN 
./RITE 

IF KOJ THEN 
BEGIN 

. DIET := CP ADS ( K JR [ J ] -RJR [ X ] , RM I [ J ’ - RP I [ X] ) ; 

• JNCERT := ABSERR [K] ; 

CEFT := D I 3T-UMCERT ; 

IF CEET>= ADSERR [ J ] THEN 
B l . G I I 

P := ?*GERT; 

1 1 := il-l; 

END; 

EuD ; 

END; 

OLDER- := ABSERR [J]; 

ADSERR [ J ] := RELERR [ J ] / P ; 

IF I>1 THEM 

ADSERR [ J j := EX? ((]. 0/M) *LM (ADSERR [J] )) ; 

T-7- r T 1 /^r on m j — m 

X — tiivwLiAi'- [J J \ O U U J 1 X\ i 

— i > r t n . 

— i iA vj i-i f 



END; 

UNTI L SHRUNr. = FALSE ; 

SC HALF; 

; (* TEND 2 *) 

PROCEDURE RPC-R (RNDEG: INTEGER; COEF : L ISTRPl ; VAR FI ERR: IN: 

VAR RJI , RJR : LISTR) ; 



EG CR ; 



AAR J , K , RLC.7 , R I GH : INTEGER; 

D EC I ? : REAL; 

A A , R VFC : MATRIX; 

RS: STRING [CO] ; 

CH: CHAR; 

“SCALE: LISTR; 
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RINT: LI ST I; 

BEGIN (* RFCR *) 

RIERR := 0; 

if ccrr[i]=c THEN 
BEG IN 

• ' - r-* r- O. 

* * - L *-• 1 '' * " • “ - t 

HR I TEL'.: ( ’ LEAL IMG 
SPACEBAR; 

T? V x T> • 

^ 1 , 

END; 

IF CCEF [1] =0.C THEN 
RECIP := i.G 
ELSE 

RECIP := I . O/COEF [ 1 ] ; 

FOR K: =1 TO RNDEG BO 
BEGIN 

RA[].,K] := -CCEF [K+l] *RECI?; 

EC 

.C; 

o.il ; 

FOR E : =2 TO REBEC 00 
RA [ K , K- I ] := 1 . 0 ; 

HAL, 'E-ICE (RNDEG, REBEC, RA , RLOvI, RICH, RSCALE) ; 
ELMHES ( RNDEG , RNDEG , RLOVj , RI GH , RA , F: I NT) ; 



!0 T? T R' t ^ 



r» -\ r j 
“ l t 



eltra:: ( an 



i< Li U »4 



Hi, RA , RI NT f RVEC) 



IIQR2 ( RRDEu , AND EG / RL0.7 , RI 311 , RPv , RXR , RA T I , RVZC , RI JAR) ; 
IT RIRRROC rHEN 

TT T? t -i 



• T O T3 

. i -i iA- 

n; T nor »f f I m 
; iv i i LLu ( _< 



= j- ; 

I MP 



I U 40 CR ITERATIONS IN XPQR 1 ) ; 



5P \CECAR ; 
EXIT; 

T? . 

Lj IV u f 

END; (* RP OR *) 



BEGIN (* TTR.PQR *) 
CLEAR5CEEZE ; 

REPEAT 

CLEARIT (1) ; 

7RITE ( ' CHAT IS TH 
INTREAD (ED EG) ; 



DEGREE OF THE POLYNOMIAL? 



n ' 



UNTIL ( ( NO EG > 0 ) A HO (NDEG<= 1AXDEG) 
PDEG := NDEG+1 ; 

FOR I : =1 TO lAEDEGPl DO 
PRCCEF [ I ] := C.O; 

r* r tt t' T 4 



EAR IT (C) ; 



iIRITELU ( ' ENTER COEFFICIENTS OF POLYNOMIAL TERMS: 
JFITE ( ' 



' ) ? 



GETPOLCC F (FDEC, PI CCEF) ; 

•'"* r ^ -* pj -,c f p p " * ‘ . 

XRITXLN ( * PLEASE XAIT ' ); 



n -n r? n y 



RPQP. (MDCo, PR 
CLEARoCREEX ; 

IP T I ERR=G THEM 



, 1 a. J_, 



r** n. T m ir \ . 

f 1 .< I f i * ) t 
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BEGIN 

NRITELN; 

WRI TELN ( ' THE ARRAY OF COMPLEX ROOTS IS'); 

WRITELN; 

FOR I : =1 TO NDEG DO 
BEGIN 

.,’RITE (’['); 

•• T-T) rnn / r- r ;p f T 1 \ . 

i u ( i .JK l i j ] , 

"• vi iu l f ) r 

ELITE (TNI [ I ] ) ; 

••’KITELN( ' ] ’ ) ; 

END; 

NRITELN; 

SPACEBAR; 

REND 2 (NDEG, PRCOEF, TER, Till , TABS EP R , TRELERR , TKLUST , TKI 
NRITELN ; 

NRITELN ( 'THE ARRAY OF ABSOLUTE ERRORS IS'); 

NRITELN ; 

’7 ~ T • — 1 ^ C • inTT r\p 

n re T ' t 

uUi . t 

mT” / i r m • 

. .. 1. 1 x-. ( l ; , 

./LITE (TABS ERR [ I ] ) ; 

WRITELN (’]’); 

END ; 

i ?n 7 n^r vt . 

‘i L\ i i i-i «j 1 1 f 



(* TTRPCTa *) 

in. 
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(* VERSION 0207 *) 
MODULE OVERLAY2 2 ; 



EXTERNAL 


PROCEDURE 


GOTOXY (X, Y 


EXTERNAL 


PROCEDURE 


CLEARSCREE 


EXTERNAL 


PROCEDURE 


GPACEGAR; 



PROCEDURE : 


:ELP ; 




TYPE STRUG* 


0 = STRING [4 0]; 




VAR I UFO: / 


iRRAY [ 1 . .12] OF STRHG4C; 




I: INTEGER; 




BEGIN (* HELP *) 




INFO [1] 


= 'TO DO MATRIX YORK ONE JUST FIRST 


r/F ,? l . 


INFO [2] 


= 'IN A -MATRIX "A". MITH THIS ONE CA 


'* » . 

f 


I iFO [2 ] 


= ' CALCULATE IT J EIGENVECTORS, EIGEN 


i . 

~ *r 


INFO [4] 


= 'VALJYS, INVERSE 01 DETERMINANT. 


t 


INFO [ 6 ] 


= 'TO SOLVE A SYSTE 1 OF SIMULTANEOUS 


l . 
/ 


INFO [7] 


= 'EQUATION'S, "A" TIMES "X" = "E" , ' ; 




I MFC [S] 


= 'ONE MUST ALSO KEY IN A MATRIX "B" 


• / 


I F C r r J 


_ i i . 

/ 




INFO [ 0 ] 


— 1 ! . 
/ 




INFO [10] 


:= 'THIS 1ATRIX "3” MAY EE KEYED IN 


•> r-» r-t — ■ i 1 . 

* L * - : ' / 


INFO [11] 


: = ' THE MATRIX "A" HAS SEEN KEYED IN 


: ND ' ; 


;‘7’' r i 


:= ' NCRKED MITH. '; 




r~* f n ^ -> Q p P rr T7 * T . 




GOTOXY ( 0 , 


- C ) ; 




NRITE (’ DIRECTIONS FOR INPUT: '); 

GOTOXY (0,2) ; 




FOR I : =1 


TO 12 DC 




; JF ITELI 


:(INFC[I] ) ; 




GOTOXY (0 ,20) : 




J FAC GEAR ; 







ELD ; (* HELP *) 



MOD END. 
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(* VERSION 0236 *) 
MODULE CVERLAY2 4 ; 



(* 


EIGEHVEC 


MCDUL 


E *) 
















TYP 


; E DCMAKII 


= 1. 


.20; 


















MATRIX 




AY [DOM A I 


MI 


,00, 


!A 


I 


L 4 


i] 


OF REAL; 




LIFT I = 


A i\ R t *» 


Y [DOMAIN 


1] 


OF 


I 


N 


r i 


EG 


ER; 




LI3TR = 


ARRA 


Y [DOMAIN 


1] 


OF 


F. 


E 




r • 

— 1 f 





VAR IONDIM: EXTERNAL INTEGER; 

IOTA , TVEC : EXTERN/ L MATRIX; 
TEVR , TEVI : EXTERNAL LISTR; 



EXTERNAL PROCEDURE CLCARSCXESM; 
EXTERNAL PROCEDURE SPACEBAR; 



:pnal function ye 



^ ^ ^ ^ r t? A ' ' • 

O. O V- V Li IN , f 



EXTERNAL [5] PROCEDURE NCR2 ( NM , N , LON , I GE : INTEGER ; VAR M : MATE I X ; Y'l 

T T • r T C^D . 

, / 1 I 1J 1 IJ 1 1 

VAR 2 .'MATRIX; VAR I ERR : INTEGER) ; 



INTERNAL [6] PROCEDURE BALANCE (NM , M : INTEGER; VAR A: MATRIX ;VA! 



: INTEGER; 

VAR SCALE: LISTR) ; 



EXTERNAL [6] PROCEDURE ELMIIES ( MM , N , LC N , I GH : INTEGER ; VAR A:1ATRIX;; 

I NT : L I S T I ) ; 



EXTERNAL [6] PROCEDURE ELTRAU ( NM , N , LON , I GH : I NTEGER ; VAR A : MATRIX? 

LI ST I ; 



Z : MATRIX) ; 



EXTERNAL [5 ] PROCEDURE 3 A LB A K ( NM , N , LON , ICH : I NT EGER ; SC ALE: LISTR;;-' 

I NTDGZR ; 

VAR Z : MATRIX) ; 



EXTERNAL [2] PROCEDURE HAFCCOPY ( HNUM : INTEGER) ; 
PROCEDURE TTRNAA ; 



VAR QUITFLAG: BOOLEAN; 

I , J , PD EG , ND IM : INTEGER; 

CHOICE, SELECT I ON: CHAR; 
TCOEF : LISTR; 

TNI , T'.JR : LISTR; 

T I ERR , MDEG : INTEGER; 

C M : C i i A R ; 

SREDL : STRING; 

P4 : INTEGER; 
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TA: MATRIX; 



PROCEDURE RNAA ( RU'D IM , RN : INTEGER; VAR RA : MATRIX ; VAP REVR , R 

VAR RVEC: MATRIX; VAR RI ERR : INTEGER) ; 



:VI : LIST 



STRICT; 

RECALL: LISTR; 

rjvrn. f T C i ” 1 T • 

i .i * I • i 1 <j i 1 ^ 

I , J , RLOiv , RICH : IUTEGER ; 

CH: CHAR; 



PROCEDURE EERMG52 ; 

BEGIN (* ERRMS32 *) 

NRITELN ( ' MORE THAN ICO QR ITERATIONS NEEDED '); 
NRITEL'I ( ' FOR SOME EIGENVALUE IN RNAA . ' ) ; 
SPACEBAR; 

EXIT; 

END; (* ERRMBS2 *) 

BEGIN (* RNAA *) 

BALANCE (END I Mi, KIT , RA , RLCU , HIGH , RSCALE) ; 

ELMHES (RNDIM, RN, R.LOw , RICH , RA, RINT) ; 

ELTRAN ( END I M , RN , RLOV , R I Cl: , RA , R I NT , RVEC ) ; 

HQR2 ( RNDIM , RM , RLCN 7 , R I GH , RA , REVR , REV I , RVEC , R I ERR) ; 

ie (rierroo) the:-; 

ERR 1CS2 ; 

BALD AX ( RNDIM , RN , RLCN, RICH , RSCALE , RN , RVEC) ; 

END; (* RNAA *) 

PROCEDURE PRVECTORS; 

VAR I: INTEGER; 



BEGIN {* PRVECTC RS *) 



- 7 

REPEAT 



:vi [ J] < > 0 . 



:rzn 



BEGIN 

NRITELN (’ VECTOR ',J,' ’HAS COMPONENTS '); 

FOR I : =1 TC NC IM DC 
BEGIN 

..’RITE ( 'VECTOR ', J ,'[', I ’ , TVEC [I , J] ) ; 
NRITELN ( ' , ' , TVEC [ I , J + 1 ],']') ; 

END; 

V.RITZLN; 

SPACEBAR; 

JRITELN ; 

NRITELN (' VECTOR ',J+1,' HAS COMPONENTS '); 

FOR I : = I TO MEIN DO 
BEGIN 



NRITE ( 'VECTOR' , J+i, ' [ ' , I , 
NRITELN ( ' , ' ,-TVEC [ I , J + l] , ' ] ' ) 
END ; 



= , TVER [I, J] ) ; 
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WRITELN; 

J : = J + 2 ; 

END 

ELSE 

T 4 L? T ■ t 

1 Li 

JRITZLH ( ' VECTOR ' , J, ' HAS COMPONENTS ') 
FOR I : =1 TO EDIM DC 
BEGIN 

JFITE ( ' VECTOR' , J I , TVS 

JR ITELE ( ' , 0 . 00 0 0 0 0 0 0 0 00 C 0 0 E + C 0 0 j ' ) ; 

TP \1 T“\ . 

J := J+l; 

EED ; 

SPACED A R ) 

CRITELN; 

UNTIL J>NDi:i; 

:e:d ; (* PR VECTORS *) 



EC I *J ( * TTFCIAA *) 

CLFADSCF EEC ; 

EDIA : = ICFDI/i ; 

FOR I : =1 TC ::dim DC 
FOR J : =1 TO ED IM DC 

TA [ I , J ] := I CTA [ I , J ] ; 

r f T? ^ . 

L i J i, i' - \ i\ L-i i—t L -J f 

J'RITE ( ' PLEASE AMT . . . ' ) ; 

REA A ( CD I M , ED I M , TA , T EVR , T EV I , T VEC , T I ERR ) ; 



( EITELN ( ' THE ARRAY CF 



JOMFLEX EIGENVALUES I! 



JRITEL”; 

FOR I : =1 TO NDIM DC 
REGIE 

J RITE ( ' r ' ) ; 

JF.ITE ( TEVF [ I ] ) ; 

WRITE 

WRITE (TEVI [ I ] ) ; 

JRITELE ( ' 7 • ) ; 

n \f p # 

Zu l>! iZ f 

JRITELE; 

SPACEBAR ; 

JRITELE ; 

AT I TELE ( ' THE COMPLEX EIGENVECTORS ARE ' 
JRITELE; 

PRVECTCR3 ; 

JRITELE; 

JRITE ( ' DO YOU WANT A HARD COPY? Y/E '); 
IF YES THEE 
HARDCOPY ( 1) ; 

ED; (* TTREAA *) 

AD DUE . 
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